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The  neutron  wave  technique  has  been  used  to  study  a  heterogeneous 
system  consisting  of  a  natural  uranium  rod  surrounded  by  a  moderator 
of  heavy  water.   By  considering  only  a  simple  unit  cell  configuration, 
the  theoretical  derivations  are  simplified  and  one  can  perform  a 
relatively  simple  and  inexpensive  experiment.   The  experimental 
technique  applied  permits  the  measurement  of  several  heterogeneous 
parameters  and  can  ba  used  to  study  ocher  fuel-moderator  systems. 

The  Feinberg-Galanin  Constant,  7,  was  measured.   This  m.easurem.ent 
was  accomplished  by  first  obtaining  the  experimental  dispersion  law 
for  the  unit  cell.   Next,  the  theoretical  dispersion  law  is  computed 
for  various  values  of  7.   The  theoretical  dispersion  law  best  fitting, 
the  experimental  one  determines  tV.e   value  of  7. 

The  theoretical  ~cdel  used  is  the  Age-Diffusion  Model.   3y  using 
a  series  of  trans form^ations  and  an  expansion  in  transverse  ei^en- 
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functions,  the  two  coupled  equations  are  solved  for  the  thermal  flux. 
Inversion  of  the  solution  yields  the  theoretical  dispersion  law.   This 
dispersion  law  is  interesting  in  that  a  critical  frequency  is  pre  - 
dieted  at  v;hich  the  dispersion  law  for  the  unit  cell  is  equal  to  that 
of  the  pure  moderator.   This  intersection  allows  one  to  measure  two 
other  heterogeneous  parameters  of  the  group  L  ,  t  ,  ,  t^  and  p. 

The  thermal  neutron  pulse,  as  a  function  of  axial  position,  was 
numerically  transformed  from  the  time  to  the  frequency  domain.   The 
frequency  dependent  attenuation  constants  and  phase  shifts  were  then 
obtained  and  the  experimental  dispersion  laws  v;ere  determined. 

The  real  and  imaginary  parts  of  the  square  of  the  complex  inverse 
relaxation  length  v?ere  expanded  in  even  and  odd  pov;er3  of  the  radian 
frequency.   From  the  expansion  coefficients  7  was  determined.   However, 
it  was  found  that  more  accurate  results  were  obtainable  when  7  was  de- 
termined from  the  dispersion  law  fit.   The  latter  analysis  yielded 
the  result: 

7  =  0.25  +  .02  era. 
The  expansion  coefficient  analysis  yielded  the  following  results: 

7(P^^)  =  0.29 


and 


rCPp  =  o.?i  . 


The   dispersion  law  intersection  was  found  to  occur  at 

Q:  =  -0.0552  +  .0003  cm" 

I  =  -0.0335  +  .0034  cm"  . 
This  corresponded  to  frequencies  of  I5S  cos  in  the  unit  cell  and 
129  cps  in  the  pure  mode.rator  dispersion  laws.   Evaluation  of  p  and 
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T  .  results  in  unrealistic  values  for  their  quantities  indicating 
th  ^ 

that  a  more  accurate  model  than  Age-Diffusion  Theory  is  required 
for  this  analysis.   Furthermore  it  is  necessary  to  eliminate 
variations  in  heavy  water  purity  to  improve  tha  accuracy  of  the 
experimental  data. 

The  analyses  of  this  work  have  been  conducted  primarily  in 

the  p  plane.   It  appears  that  more  accurate  and  sensitive  results 

2 
are  obtainable  when  the  analyses  are  conducted  in  the  p  plane. 

Data  to  aid  this  type  of  future  studies  have  been  included. 
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CHAPTER  I 
INTRODUCTION 

Background 

The  Department  of  Nuclear  Engineering  Sciences  at  the  University 
of  Florida  has  been  engaged  in  an  extensive  research  program  in  the 
field  of  nevitron  wave  propagation  since  1960.   The  original  studies 
dealt  vjith  the  theoretical  (1)  and  experimental  (2,  16)  determina- 
tion of  the  diffusion  and  therraaliza tion  parameters  of  homogeneous 
moderating  media.   Successful  results  prompted  studies  of  two  region 
rioderating  media  (5),  moderating  madia  with  localized  absorbers  (4), 
a  "homogenized"  subcritical  assembly  (5,  6)  and  a  reflected  subcritical 
assembly  (7,  3,  11) .   The  continued  success  of  the  neutron  wave 
technique  has  been  demonstrated  in  recent  conferences  held  at  Karlsruhe 
(9),  Gainesville  (10),  San  Diego  (11)  and  Ann  Arbor  (12).   Since  oast 
endeavors  have  bean  successful  it  is  of  interest  to  investigate  the 
application  of  this  technique  to  other  syste^.s. 

In  this  work  a  heterogeneous  moderator-fuel  system  will  be  investi- 
gated using  the  neutron  wave  propagation  technique.   Cain  (15)  studied 
the  heterogeneous  problem  theoretically  and  com.puted  the  spatial  and 
frequency  dependent  neutron  fluxes  in  various  configurations.   One  of 
the  configurations  consisted  of  two  natural  uranium  fuel  rods  embedded 
in  a  heavy  water  moderating  medium.   Calculation  of  the  flu:<  per -.urba  tion 


2 
showed  no  significant  change  due  to  the  fuel  rods.   Experimental 
measurements  were  predicted  to  be  indeterminate  unless  enriched  fuel 
was  used.   Measurement  of  quantities  other  than  the  flux  would  have 
to  be  made  to  obtain  meaningful  results  with  natural  uranium  fuel  rods. 

Corno  ilk)    also  studied  heterogeneous  systems  theoretically  and 
proposed  an  experiment  to  study  the  heterogeneous  characteristics  of 
a  cylindrical  fuel  rod  embedded  in  a  moderating  m.ediun.   The  experi- 
ment involved  the  use  of  a  highly  enriched  driving  shell  surrounding 
the  fuel  rod.   By  measuring  the  exponential  attenuation  of  the 
fundamental  mode  for  different  driving  shells,  7  and  1^  can  be  determined. 
7  is  the  Feinberg-Galanin  Constant  (15)  and  t^  is  the  nunber  of  neutrons 
produced  per  thermal  neutron  absorption  in  the  fuel.   The  technical 
problems  involved  in  designing  and  using  such  a  driving  shell  maks  a 
simpler  experitr.ental  configuration  desirable. 

Booth,  Perez  and  Hartley  (16)  have  developed  a  procedure  for 

neutron  wave  propagation  exper  ir.ents .   This  procedure  is  a  general  one, 

applicable  to  aln-.ost  any  system.   The  experimental  quantity  measured 

is  the  cor.plex  inverse  relaxation  length,  p,  as  a  function  of  frequency 

where 

p  =  2  +  i|, 

at  is  the  amplitude  atten'jation  per  unit  length  and  |  is  the  phase 
shift  per  unit  length.   By  measuring  p  it  is  then  possible  to  obtain 
the  dispersion  law,  i.e.,  :  versus  c:  with  frequency  as  the  para-eter. 
Analysis  of  the  complex  inverse  relaxation  length  allows  the  experi- 
mental de  zeriT'.ina  tioa  of  the  Feinberg-Galanin  Constant  for  the  fuel  rod 
in  the  syste-.   A  ccmaarison  bef-'een  Che  dispersion  lav  predicted  by 
thc-ory  and  the  experir.enta]  dispersion  law  can  also  be  rT-.ade.   This 
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method  of  analysis  will  be  employed  to  study  the  heterogeneous  system. 

Scope  and  Qbiectives 

The  neutron  wave  technique  has  been  used  successfully  to  study 
several  systems.   It  has  not  been  used  experimentally  to  study 
heterogeneous  moderator-fuel  systems  and  the  question  of  its  practi- 
cality in  such  cases  remains  unresolved.   This  study  will  investigate 
the  use  of  the  technique  on  a  system  other  than  those  already  con- 
sidered.  The  scope  of  the  neutron  wave  technique  is  to  be  extended 
experimentally  to  include  a  method  for  investigating  the  heterogeneous 
moderator-fuel  system.   The  system  is  to  be  as  simple  as  possible. 
Heterogeneous  parameters  and  their  sensitivity  to  the  neutron  wave 
parameters  will  be  determined.   In  particular,  the  Feinberg-Galanin 
Constant,  7,  will  be  measured.   One  natural  uranium  fuel  rod  embedded 
in  a  moderating  medium  of  heavy  water,  as  shown  in  Figure  1.,  will  be 
used  to  illustrate  the  method. 

Two  other  heterogeneous  parameters  of  the  group  t^,  p,  -r  ,  and  L 
can  also  be  measured.   Their  measurements  are  made  possible  by  the 
intersection  of  th.e  dispersion  laws  of  the  moderator  only  and  the 
modera tor- fuel  systems. 

A  review  of  the  literature  failed  to  show  the  existence  of  an 
established  experimental  method  to  measure  7  directly.   7  is  usually 
computed  but  complex  designs  oftp.n  call  for  approximations  which  lead 
to  uncertainties  in  the  results.   The  method  presented  here  allows  the 
measurerrent  of  7  to  be  made  in  a  si-^.ple  and  inexper^sive  way  using 
one  fuel  rod. 


FIGURE 


A  sche:-u\tic  vie'W  of  the 

HETEROGENEOUS  SYSTEM 


The  effect  on  7  when  a  lattice  is  considered  needs  to  be  examined. 
Direct  extrapolation  from  the  one  fuel  rod  case  to  that  of  a  complex 
lattice  must  be  demonstrated  experimentally.   In  the  absence  of  such 
data  an  extrapolation  to  the  simple  lattice  configuration  used  by 
Dunlap  (5)  will  be  investigated.   If  lattice  effects  are  significant 
the  experimental  method  must  be  modified  to  make  the  measurement  of 
7  a  valid  one. 

The  potential  of  the  method  will  be  investigated  and  its  exten- 
sion to  determine  other  heterogeneous  effects,  such  as  interaction 
and  shadowing  effects,  will  be  considered.   Versatility  of  the  method 
v;ill  also  be  discussed. 


CHAPTER  II 
THEORY 

Introduction 

The  theoretical  model  used  in  this  analysis  is  the  Age- 
Diffusion  Model.   Although  this  is  a  sinplified  approach  it  is 
sufficient  to  exhibit  the  characteristics  to  be  expected  from  a 
more  refined  analysis.   To  be  more  accurate,  thermalization 
effects  should  be  considered.   However,  Dunlap  (6)  found  that  ex- 
periraental  results  for  the  suLcritical  assembly  agreed,  within 
several  percent,  with  Age-Diffusion  Theory  for  frequencies  up  to 
330  cps.   .As  will  be  seen  in  the  development  of  this  chapter, 
the  theoretical  basis  for  the  experimental  method  is  established 
by  considering  trequencies  less  than  330  CuS.   Thereicre  the 
Age-Dittusion  Model  will  be  used. 

The  purpose  ot  the  theoretical  model,  in  this  case,  is  to 
predict  values  for  param.eters  which  are  to  be  measured  experi- 
mentally.  .Age-Diffusion  Theory,  since  it  does  ignore  thermaliza- 
tion effects,  can  not  be  expected  to  predict  accurate  values. 
However,  these  predictions  are  sufficient  to  verify  the  practicality 
of  the  experinant  and  do  offer  a  comparison  with  the  experimental 
results.   Better  comparisons  will  be  obtained  with  more  accurate 
theoretical  models. 


An  important  nurr.ber  associated  v;ith  heterogeneous  calculations  is 
the  Feinberg-Galanin  Constant,  7  (15).   This  constant  is  defined  as 
the  ratio  of  the  net  current  cf  neutrons  into  a  fuel  rod  divided  by  the 
neutron  flux  on  the  surface  of  the  rod.   The  value  of  7  is  determined 
by  tv;o  methods.   The  first  method  is  to  fit  the  experimental  disper- 
sion law,  using  7  as  the  fitting  parameter.   The  second  method  is 

2 
based  on  expanding  p   in  the  dispersion  law  in  a  power  series  of  ioo 

and  equating  coefficients  of  like  pov;ers.   The  first  two  expansion 

coefficients  will  be  related  to  7.   The  results  of  the  two  methods  will 

be  con^.pared. 

In  the  theoretical  development  of  this  chapter  it  will  be  shown 

that  the  dispersion  laws  for  the  moderator  only  and  the  moderator-fuel 

rod  intersect.   Tht^  frequency  of  the  moderator-fuel  rod  dispersion  law 

at  which  the  intersection  occurs  is  designated  the  critical  frequency. 

CO  .   Exoerimsntallv  this  is  important.   By  setting  the  real  and 
c 

imaginary  parts  of  the  moderator-fuel  rod  dispersion  law  equal  to  those 
of  the  moderator  dispersion  law  and  substituting  the  value  of  the 
measured  critical  frequency  and  Zt   and  |  measured  at  the  critical  fre- 
quency, two  heterogeneous  parameters  can  be  deceruined.   Further  de- 
tail is  presented  later  in  this  chapter. 

The  Theoretical  Model 

Two  equations  make  up  the  Age-Diffusion  Model.   The  time  dependent 
fast  flux,  6  (k,  u,  t)  ,  in  tert-.s  of  lethargy,  u,  is  described  by  the 

equation 

-D(u)  V""  «(x,  u,  t)  +  £^(u)  ^(x,  u,  t)  - 

[2.1] 
d^(x.  u.  t)  ^     S(x,  u,  t)  _  Sq(x.  u,  t) 


v(u)      dt  ^ 
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and  the  associated  equation  for  the  thermal  flux,  >>   (x>t)  is 


^th  2   th,   ,    „th   th, 

D   V  4   (x,t)  +  I,   0   (x,t)  = 


:^  th 
1 d<> 


[2.2] 


^th 


^^    +  q(x,a^^,t)  -  A(x,t) 


where 

D  =  the  diffusion  constant  [cm] 

Z  =  the  macroscopic  absorption  cross  section  [cm  ] 

,   .    r       -li 

V  =  velocity  [cw  sec   J 

q  =  slowing  down  density  [neutrons  cm   sec   ] 

S  =  source  of  neutrons  due  to  the  fuel  rod  [neutrons  cm   sec   ] 

A  =  absorption  of  thermal  neutrons  in  the  fuel  rod  [neutrons  cm  "^sec   ] 

Since  cylindrical  geometry  will  be  considered, 
X  =  x(r,  z) 

in  both  equations. 

Solution  of  Equation  For  Fast  Neutrons     The  time  dependence  in 
equation  [2.1]  is  removed  with  a  Fourier  Transformation  with  respect 
to  t.   Since  the  flux  results  from  the  pulsed  neutron  source  which  is 
assumed  to  stare  at  zero  tir-.a 

*(x,u,t)  =0;   t  <  0. 
Under  this  condition  the  Fourier  Transform  of  equation  [2.1]  is 

2 

-   D(u)   V     <>(x,u,'jj)    +      L   (u")    i>(x,u,05)      = 

^     '  [2.3] 

where  o)   is    the    transform  parameter. 


The  source  term,  S(x,u.a3),  represents  the  fission  neutrons  produced 
by  thermal  neutron  absorption  in  the  fuel  rod.   It  is  assumed  that 
these  source  neutrons  appear  with  lethargy  zero,  at  the  surface  of 
the  fuel  rod.   The  explicit  form  of  the  source  term  is  then 

c/            N     =^  /  N  ^  (o)     5(r  -  r  )   th, 
S(x,u,a3)  =  5(u)   s''   T^  7  _J ol  *   (x,u)) 

Z  (o)      2  ir  r 


where 


r  =  radius  of  the  fuel  rod  [cml 
o  '■       ' 

s  -  the  probability  of  a  neutron  entering  the  slowing 

Z  (o)  down  process 

Tj  =  the  neutron  yeild  per  thermal  neutron  absorption 

7  =  net  thermal  neutron  current  into  the  rod        [ cm] 
thermal  neutron  flux  at  the  surface  of  the  rod 

5  =  the  Dirac  Delta  Function. 


It  is  convenient  to  write  equation  [2.3]  in  terras  of  the  slewing  down 
density,  q(x,u,cD).   The  infinite  medium  relation, 

q(x,u,a^)  =  :  Z^(u)  5(x,u,a^), 

is  assumed  since  tr-easurements  were  made  in  heavy  water  which  is  a 
weak  absorber  of  neutrons.   Substituting  for  i(x,u,a)  and  S(x,u,a)) 
equation  [2.3]  becomes 

{  -  D(u)  V^  ^  Z  (u)  +  ^  }  ^-^^-   = 
a      v(u)     I   '^.(u) 

So(x,u,a.^     ,,  ,  ^s^°>       ^^'^  -  ^o^   th,   , 
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A  change  of  variables  is  now  made  from  lethargy,  u,  to  age,  t,  using 
the  definition 


du 


Jli}± 


^(t1 


and   the   relations 


du" 


J>iT). 


i    ^t^T) 


5(u)   = 


D(t) 


4    ^t(T) 


6(t) 


equation  [2.1i.]    becomes 

(    -  V^  +   PCt.o.)}    q(x,T.a))    =   -      M|^T^ 
_,   .       E   (o)  ^      5(r-r   ) 


^Jo) 


2  rt   r 


th, 
0      (x,a)) 


[2.5] 


where 


PCt.gj)    = 


D(t) 


By  multiplying   equation  [2.5]    by    the    integrating   factor 


o 


p(t,uj) 


icoL 
p  e  s 


where  L   is  the  slowing  down  time  and  o  is  the  resonance  escape 
probability  and  recognizing  that 


^T      P(t,^)         PCt.'JJ)       3t  P(T,a3)  ' 
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p       ?^  ^Jo)      •<\y         5(r  -  r  ) 

{-V    .  ^-  }  e(x..,a.)  =    5(v)  -^  -(^;^   -YT7-  *    (-'^)    t2.6] 


where  , 


Equation   [2.6]    ccmpares   with  equation  [8.2.9a]    of  Beckurts   and 
WirtE   (17). 

The   next   step   is    to   remove    the   age   dependence    in  equation   [2.6] 
using  a  Laplace   Transformation  with  respect    to   t-      The   boundary 

conditions 

e(x, -7=0,00)    =   0 

and  p(T=o,aj)    =   1 

apply  since  no  neutrons  slow  down  past  age  zero  and  the  resonance, 
escape  probability  at  age  zero  is  one.   The  former  condition  holds  be- 
cause the  source  neutrons  are  taken  into  account  by  the  source  term 
in  equation  [2.6].   The  Laplace  Transformation  of  equation  [2.6]  is 

^2  ,     .      ,     .     E  (o)        5(r  -  r  )^th.   .   '  „  „, 
V  G(x,s,gd)  -  se(x,3,co)  =-  s^    T]  7 o'<t>   (x.o))   L2.7] 

E  (o)~         2  n  r 


where    s    is    the    transform  para-.eter. 

At    this   point    the   dependent  variables,    9(x,s,cd)    and  (S      (x,a5)  , 

are   expanded   in    terms   of    the   radial   eigenf unctions ,    J    (Br),    given 
'^  °  o      n   o 

by   the    solution   of    the  Kelinholtz   equation, 

(V^  +  B^)    J    (B   r)    =   0. 
r  n        on 


The  orthogonality  condition  is 
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R 
o   o   n    o   p  p  -"np 

where  N   is  the  normalization  constant  and  5   is  the  Kronecker  delta, 
R,  the  extrapolated  boundary,  is  such  that 

J^(B^R)  =  0. 

The  expansions  are 

G(x,s,a))  =  E  e  (z,s,cd)  J  (B  r) 
n   n         on 


and 


.'\x.c.)  =  L  ^Jz,..)    J^(B^r) 


By  noting  that 

r  J  (B  r)  =  -  B"  J  (B  r), 
r   o   n        n   o   n   ' 


V^  J  (B  r)  =   "^ 


multiplying  by 

2  rt  r  J  (B^r). 
o   p 

integrating   and  applying    the   orthogonality  condition,    equation   [2.?] 
is    reduced   to    the    following   equation   for    the   expansion  coefficients: 


2 

^  e  (z.s.cu)        .0 

°  -    (B      +   s)    9    (2,s,c^)      = 

V    2  P  P 

dz 


[2.8] 


fs^       IZ     L     /^    (z.o.)    J    (P    r   )    J    (B   r   ). 
Z    (o)         N  ""        ^  o      n   o        o      p    o' 

t  P 


The  last  step  in  the  solution  to  the  fast  flux  equattor.  is  to  rs- 
move  the  spatial  dependency  by  a  second  Laplace  Transformation,  this 
time  with  respect  to  2.   The  Laplace  Transformation  of  equation  [2.6] 


is 
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2  2 

P  e  (p,S,co)  -(B   +  s)  G  (p.s.cd)  = 


^s^^^  "^^        L^^^    (p/x)  J  (B  r  )  J  (B  r  ) 
Ti^  F~   ^   "        o  n  o   o  p  o 
t^  ^   p 


where  p  is  the  transform  parameter  and  the  boundary  conditions  are 
taken  to  be 


d9  (z=o,s,a)) 


P 

e  (z-o,s,a))  =   — ^, .  =  0 


which  neglects  the  gradient  at  z-o  because  the  dispersion  lavr  is 

independent  of  its  presence.   The  solution  for  the  fast  flux  expansion 

coefficients  is 

[2.9] 

%<p'='-)  -    — —2 — rt^j  -¥-  ^\  (p'-)  -^oVo^  -^o^^V- 

^  s  +  (B_^  -  p  )   t^  '     p 


In  order  that  the  solution  may  be  of  use  in  solving  for  the  therraal 
flux,  it  is  necessary  to  transform  from  the  s  domain  back  to  the  t 
domain.   The  inverse  transformation  of  equation  [2.9]  is 

2    2 

P     r,  (o)   n  7    ^.  th, 


eJp,T,co)  =   e  IaI::^   11  Z  o:'-(p,a))  J„(3^rJ  J^(3^rJ 


p  IToT    ir~  "  "  ■       °  ^  °    °  p  ° 


n,p  =  1,2,  .  .  .  k. 

Solution  of  Equation  for  Thermal  Neutrons    It  is  also  necessary  to 
solve  the  thermal  flux  equation,  equation  [2.2].   The  procedure  is 
similar  to  that  used  to  solve  the  fast  flux  equation.   The  order  o£ 
transformation  is  cn:inged  so  that  the  boundary  conditions  may  be  more 
readily  applied.   Before  solving  equation  [2,2]  it  is  noted  that 
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* /   N     6(r  -  r  )     th, 
A(x,t)  -  7  _J ol   4.   (x,t) 

2  It  r 
which  represents  the  absorption  of  neutrons  in  the  fuel  rod  as  losses 
occurring  at  the  surface  of  the  fuel  rod.   By  representing  the  absorp- 
tions in  this  manner  and  the  source  term  in  the  manner  described  in 
the  previous  section,  the  problem  is  reduced  to  a  one  region  problem 
with  a  source  and  absorber  of  neutrons. 

Equation  [2.2]  is  Laplace  Transformed  with  respect  to  the  z 
variable: 


[2.10] 


,  _th  2   ,,th  2   _th-,   th.     .    ^th   th. 

(-D  p  -  D   V  +  Z   )  0   (r,p,t)  +  D  p  0   (r,z=o,t) 

.  n*^^  -^  (r,z-o,t)  =  J^  bo^     (r,p,t)  .  „,  ^  ^   .x 

th 

^th.     ,  S(r  -  r  ) 
-  70   (r,p,t)  __\ ol 

2  jt  r 


where  p  is  the  transform  parameter.   In  this  case  neither  the  thermal 
flux  nor  iti  derivative  is  zero  at  the  origin.   These  terms  are  ex- 
pressed as 

b^   (r.z^o.t)  =  foill^^  =  S  (t)  J  (B^r) 
dz  2       °     °   ^ 

and 

o'^'^Cr.z^o.t)  =  r  c  (t)  J  (B  r) 
nn     on- 

Making  these  substitutions,  using  the  eigenfunc tion  expansions  and 
invoking  the  orthogonality  condition  as  was  done  in  the  solution  of 
the  equation  for  the  fast  neutrons  equation  [2.10]  becomes 


(-dV  *  d'^B^  .  .f  )/\p.t)  .  ,i\i,,    .  B'"s^(t)5p^ 
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do 


th 


^h    "'^ 


(p,t)  +  qp(p,T^^,t) 


-2_  E  0^^\p,t)  J  (B  r  )  J  (B  r  ) 
jj     n  n        o  n  o   o  p  o 

P 


where 


q(r,p,Tj.^.t) 


^  ci  (p,T^,  .t)  J„(B  r) 
n   n     th     on 


The  time  dependence  is  transformed  out  by  a  Fourier  Transforma- 
tion, again  using  to  as  the  transform  parameter.   The  result  is 


{-p2  +  B^  + 


_th  ^    . 

Ith)    *^^    (p.a^)    +   p  Cp(u:)    +  S^Cjq)    5^^ 


=     ^-J—  Z   >^^'(p.J^)    J    (3    r   )    J    (B   r   )    +   p(t.,  .a:)    &    (c,t.,  ,'i3) 
th^      nn  ono        opo  th    p  tn 


.th 


where 


qp(p,T^.^.a»    =   p(T^.^_,a))    e^{p,^^^,-^) 
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Combination   of  Fast  and  Tberml  Neutron  Solutions        The    solution   for 
the    thermal    flux   expansion  coefficients    is 


L/r2         2         _1  iio    -1     th,        ^ 

1%    ^^    ^2-     -~'^I%    ^P''"^ 


th 


o 


P  [2.11] 


X  g  C^P'^^^-^o^Vo^-^o^^'^o^^   P  %^^>    ^  ^^^^)    ^1 


n,p  =   1,2,    .    .    .    k 

where    the    solution  obtained   for    the    fast    flux   expansion  coefficients 
has   been   substituted, 

_a =       1 

th 


and 


th  o 


Since  the  moderator  to  fuel  ratio  considered  is  large,  the  probability 
of  a  source  neutron  being  absorbed  on  its  first  collision  is  very 
small.   Therefore,  the  approximation 


is  made 
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Numerical  Solution 

Two  cases  of  tlie  solution  will  now  be  considered,  n,p  =  1  and 
n,p  =  1  and  2.   The  former  case  represents  the  solution  for  the 
fundamental  spatial  mode  while  the  latter  illustrates  the  case  which 
includes  the  first  spatial  harmonic  as  well  as  the  fundamental  solu- 
tions.  The  latter  was  used  in  all  the  computations.   The  importance 
of  the  first  harmonic  will  be  examined. 

In  both  cases  an  inverse  Laplace  Transformation  is  required  in 
order  to  obtain  the  solution  in  z  space.   To  simplify  the  inverse 
transformation  the  poles  of  the  transformed  solution  can  be  computed. 
The  solution  in  z  space  is  then  of  the  form 

/Nz.cn)  =  2G.('^)e^i^^"^^ 
P  1-   1 

where   G .  (ca)    is   a   constant,    depending   on    the   pole    p .  (co)  .      The    final 
solution   in    ter"s   of    z  will   not  be   computed.      Instead,    the   poles   of 
the    transformed   solution  are   computed,    yielding    the    dispersion   law. 
For   n,p   =    1,    p      is   computed   from   the   equation 


1      ^'l  '^^    ".2  -.    th/ 


L\  D 

th  o 

-(B        -    p   )t  ,     -    WjL  p 

+     ■    7     .    /nPe       ■  ^   -    ll    J,    (3,0    =   0                     [2.12] 
D      N 

For   n,p   =   1  and   5,    p  is   computed  by   setting    the    determinant   of   a 
2X2  matrix   ecual    to   zero.      The   elements   of    che   matrix   are 


18 


2  2 

A        =-B        +   p     -     _1_  -  _M^ 

th  o 

U'\    '^  J      o        1   o 


2  2 

.         -(B,       -   P  )t  .     -    ix.L 

^12  =  -2. inPe  -    1  ,►    J    (B   r  )J    (E   r  ). 

■••^  ,.|^  1^  J        o      1    o      o      2   o 


^2 


2  2 

"(B        -   p   )t  ,     -    iaJL 

2 


2  2 

A        =   -B        +   p      -  L_     ■  .  i^  ,  , 

th  o 


"(B        -   p   )t        -    iii3L 
■^  7  Tnpe        2  th  s      _    ^.^  J  2(B        ) 


d'\ 


J       o    ^   2   o' 


In  each  case   solutions   ware   obtained  using    the  FERVE2  Code    (See 
Appendix   D) . 

Results   of    the   computation   for, selected   frequencies   are   given   in 
Table    1.      The    fundanental   eigenvalues   changed   in    the    fourth   signifi- 
cant  digit    froiT.   the   one    tern   (n,p   =   1)    to    the    two    term   (n,p   =   1   and  2) 
cases.      Computations   were   also   -jda    for   a    ten    term  expansion    (21)    and   th': 
fundamental   eigenvalues   differed  by   less    than   3.5%   frora   those   of    the 
one    term  case.      In   the    two    terai  case    the   eigenvalues   of    the    first 
harmonic   mode   are   of    the   sarae   order   of  magnitude   as    the    fundamental 
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eigenvalues.  However,  the  experimental  source  was  found  to  contain  a 
low  percentage  of  first  harmonic  contamination.   These  results  are 
included  in  Table  1.   In  this  case  only  the  fundamental  mode  is 
significant. 


7  Analysis 

Using  equation  [2.12]  it  is  possible  to  analyze  7.  7  is  varied 
parametrically  and  the  resulting  dispersion  law  is  computed  and  com- 
pared v;ith  the  experimental  one.   The  best  comparison  represents  the 

measured  value  of  7  within  the  limitations  of  the  theoretical  model. 

2 
If  p  in  equation  [2.12]  is  expanded  in  a  power  series  of  the 


form 


p^  =  E  P  (lj)" 

'^  n  n 


and  coefficients   of    like   powers   are   equated,    still  another  method  of 
determining  7    is  available.      The    first    two  coefficients  of    the  ex- 
pansion are   related   to  7  as    follows 


_i_    +    7    X(l     -     T]    ?) 

2        L^ 
Pn    =   S,      +  \h_ [2.13] 


0  1 

1   +   q    p    7   X    -; 


th 


1.   These  results  were  obtained  frci:"!  the  expansion  coefficients  of 
the  transverse  continuous  mode  dat?  at  z  =  i+0  ca.   The  UFNLLS  Code  (20) 
was  used  to  fit  the  data  to  the  fundarr.ental  and  first  harnonic  modes. 
Contamination  at  z  =  50  cm,  the  first  data  point  analyzed,  was 
determined  by  evaluating  the  exponential  attenuation  of  each  mode  and 
computing  the  percent  contribations.   The  percentages  are  based  on 
maximum  amplitudes  and  represent  maximum  contamination. 
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and  __^,x(,pL^(l  -B^  ^,^^^,'^,^)) 


vhere 

X  =  _o_i_e  . 

These  expressions  were  obtained  by  separating  the  exponential  term  into 
frequency  dependent  and  independent  terms,  expressing  the  result  as  a 
product  of  exponentials  and  expanding  each  product  in  a  MacLaurin  Series. 
The  sensitivity  of  f   and  P   to  y   was  found  to  be 

•^-v  3-5  X  10 

and 

■^-  3.1  X 10  ^ 

both  P  and  P  changing  about  the  same  amount  in  the  second  significant 
digit  with  respect  to  y ,    yielding  a  257.  change  in  P  and  JX   change  in  ?  , 

Con;parison  of  the  two  y   analyses  with  the  experimental  results 
will  be  made  in  Chapter  IV. 

The  Critical  Freauancy 

To  investigate  the  critical  frequency  equation  [2.12]  is  again 

considered,  both  for  the  moderator  only  and  the  moderator- fuel  rod 

systems.   The  right  hand  sides  of  both  equations  are  zero  so  the  left 

hand  sides  can  be  ecuated.   L  ,  and  D  "  are  constants  of  the  moderator 

th 

but  D    differs  between  the  two  systems  because  v    is  different. 
o  ^ 

2 

However,  this  d5.fterence  is  small    and  they  are  assumed  to  be  equal. 

It  is  new  assumed  that  the  disoersion  laws  intersect 


th 


2.   Subsequent  analysis  shewed  the  difference  in  v    to  be  less 
than  2'"o. 
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so  that  the  two  p's  are  equal.   The  relation  is  then  obtained: 

in,  -/J   ^B,r   )     ,       -(^^   -   Pc^^^h   ■    ^c^  >  .      ,,, 

-  __c     +       010"^    Iripe  -    1  L  [2.15] 

o  I 

-    ico 
m 

o 
vhere  cu   is  the  frequency  of  the  moderator  only  dispersion  law  asso- 
ciated with  the  intersection  and 


c    c    *c 

IS    Liie   vctiut   <jf   p  at    the    intersection,      u)      is    the   critical    frequency 
of    the   moderator-fuel   rod   dispersion   law  at    the    '.ntersection.      The 
dispersion   laws   intersect  at    the   same  values  of  a  and   |  but    the 
associated   frequencies  are   different.      By  equating  real  and   imaginary 
parts   of   equation  [2.15]    the    two   equations 


n   p   cos    (Pa    ;    T^,     -  O)  L  )i 
'  c   c    th  c    s 


-     i\^-     C\'-5/)iT,h 


=       1 


and 


■(Bf-cv  -  i^^u^^    ^^  .^    ^th^ 


T1   p   sin(2Ci   ?    T  ,     -  'JJ  1-  )e  =     c  m     1 

'  c-c    ta  c    s  TV —     r- 

o  o        L   o 


are   obtained.      Computations   performed   for   many  conditions   of   inter- 
section  showed    that     cu     >  00    .      Therefore,    the    sine   and  cosine   must  be 

c  m. 

positive  for  a  solution  to  exist  because  all  the  other  quantities  in 
both  equations  are  positive.   The  angle 
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2Q!    P     T    ^     -    CD   L 

c^c    th  c   s 

must  be    in   the    first   quadrant.      If    these    two   equations   are   divided, 
one   by   the   other,    the   result   is 


a>  -  en  D     N,  , 

tan    (2a  E    T  ,    -  CD  L   )    =     -?— -B 1 •  [2.16] 

c^c  'th  c    s-^  ^th  -   2,„ 

D  7   J      (B.r   ) 

o  o        1   o 

By   determining    the   experimental    intersection  and  computing  N      and 

2 
Jq    ('^I'^o^    ^^'^    right   hand   side    of   equation   [2.16]    car.   ba   evaluated, 

.  -   th 

If  V   and  7  are  known.   The  argument  of  the  left  hand  side  can  then 

be  determined.   This  allows  either  t  ,  or  L   to  be  evaluated  in 

th     s 

terms  of  the  remaining  parameters.   Once  the  argument  is  found,  either 
of  the  tvo  original  equations  can  be  used  to  solve  for  t]  or  p  in 
terms  of  the  other.   The  experimental  results  will  be  discussed  in 
Chapter  IV. 

The  intersactiop.  of  the  dispersion  laws  is  illustrated  in 
Table  2.  and  shown  in  Figure  2.,  where  they  are  plotted  for  7=  0.0 
(pure  moderator)  and  7  =  0.6.   These  values  ware  selected  to  give 
a  qualitative  description  of  the  intersection.   Based  on  Cain's 
results  (15),  the  value  of  7  was  expected  to  be  in  this  range. 

The  sensitivity  of  cd  to  7  was  investigated  by  using  the  two 
term  expansion  of  equation  [2.12].  It  was  found  that  a  change  of 
0.1  in  7  produced  a  2  cps  change  in  cd  . 


24 


TA3LE  2. 

NUMERICAL  ILLUSTRATION  OF 
INTERSECTION  OF  DISPERSION  LAWS 
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CHAPTER  III 

EXPERII-1ENT.\L  METHOD 
Introduction 

Two  experiments  were  conducted.   The  first  was  in  a  pure 
moderating  medium  of  heavy  water.   The  second  was  in  the  same  system 
with  one  fuel  rod  inserted  as  shcv;n  in  Figure  1.   The  experimental 
dispersion  laws  for  both  systems  were  determined  using  the  pulse 
propagation  technique  (16)  and  converting  the  data  into  the  frequency 
domain  in  the  usual  way.   A  Texas  Nuclear  Corporation  Neutron 
Generator  (Model  9505)  was  used  as  the  source.   Tne  system  in  which 
the  measurements  were  T^ade  and  the  neutron  generator  which  was  used 
have  been  described  in  detail  by  Dunlap  (5,  6).   The  17  inch  graphite 
stack  between  the  light  water  thermalizing  tank  and  the  heavy  water 
tank  was  removed.   This  was  done  to  reduce  the  width  of  the  thern^al 
neutron  pulse  and  thus  increase  the  high  frequency  content. 

T^jo  data  acquisition  systems  were  used,  one  heaving  a  movable 
de tec  cor  which  was  used  for  measurements  along  the  axial  length  of 
the  assembly  and  the  other  being  the  reference  detector  which  was 
used  for  nor.r^lization  purposes.   The  reference  detector  was  moved 
to  a  position  adjacent  to  the  graphite  in  the  thermalizing  assem.bly. 
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This  was  done  to  reduce  the  effects  of  source  anisotropy.   A  brief 
sketch  of  the  set  up  is  shown  in  Figure  5.   Each  detector  system 
contains  basically  the  same  components  with  the  following  exceptions: 

(a)  The  laovable  detector  system  has  a  12  inch  He3 
detector  and  the  reference  detector  system  has 
a  6  inch  He3  detector 

(b)  The  movable  detection  system  has  additional 
components  for  data  acquisition  and  analysis 
of  the  neutron  pulse  in  time. 

The  two  detection  systems  are  shown  in  Figures  k.   and  5. 

Considerable  tiine  was  spent  in  becoming  familiar  with  each 
and  every  component.   This  step  can  not  be  overemphasized  in  im- 
portance for  it  proved  invaluable  in  the  iinmediate  identification 
of  equipment  failures.   Several  failures  occurred  and  data  acquisi- 
tion was  terminated  with  a  minimum  collection  of  faulty  data.   The 
use  of  oscilloscopes  was  likewise  important  in  the  verification  of 
proper  pulses,  monitoring  signals  and  general  trouble  shooting. 

Light  water  ccntamination  in  the  heavy  water  moderator  changed 
between  the  two  experiments.   The  first  experiment  was  conducted  at 
99.5"  purity  while  the  second  one  was  conducted  at  99. Oo  purity. 
The  decrease  in  purity  was  ca-jssd  by  the  several  transfers  of  the 
moderator  that  took  place  between  the  oerforniance  of  the  two 
experim.ents. 


1.   The  angular  distribution  of  neutrons  was  found  to  be  de- 
pendent on  the  amplitude  of  the  target  current.   Experinieatai  data 
shoved  that  the  variation  of  anisotropy  for  the  range  of  target 
current  amplitudes  used  was  reduced  to  less  than  1%   when  the  reference 
detector  was  positioned  adjacent  to  the  graphite. 
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FIGURE  3. 

EXPERIiMEKTAL  ASSEMBLY 
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MOVABLE  DETECTOR  COUNTING  SYSTEM 
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FIGURE    5. 
REFERENCE    DETECTOR    COUNTING    SYSTEM 
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Preliminary  Procedures 
Before  any  data  were  taken  the  following  pre-operational  checks 
and  tests  were  performed  on  the  counting  systems: 

(a)  conpatability  of  output-input  signals  between 
components 

(b)  detector  plateau  determination 

(c)  chi-squared  checks 

(d)  resolution  time  determination 

(e)  garnma  ray  discrimination  checks 

(f)  minimizing  of  internal  noise 

(g)  reverif ication  of  plateaus,  chi-squared  check 
and  resolution  time. 

Step  (g)  is  necessary  to  assure  that  changes  made  in  the  system  in 
the  preceding  steps  have  not  affected  the  counting  characteristics 
of  the  systems.  Adequate  guidance  in  the  completion  of  most  steps 
is  easily  obtained  from  the  appropriate  technical  manuals  and  any 
Nuclear  Engineering  laboratory  m.anual.  Resolution  time  determina- 
tion is  described  in  Appendix  A. 

The  last  task  to  be  accomplished  before  the  actual  experiment 
is  begun  is  the  selection  of  the  target  current  pulse  width  and  re- 
petition rate  for  the  neutron  generator.   Coincident  with  this  is  the 
selection  of  a  compatible  cViannel  width  for  t'ne  multichannel 
analyzer.   This  procedure  is  described  in  Appendix  B. 

In  the  course  of  checking  out  the  detector  systems  a  lim.iting 
count  rate  was  observed.   Increased  source  intensity  produced  no  in- 
creased scaler  count  rate  at  appioxinately  200,000  and  2^0,000 
counts  per  second  for  the  reference  and  movable  detector  amplifiers 
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respectively.   All  data  were  checked  to  insure  that  the  maximum 
count  rate  during  the  neutron  pulse  peak  did  not  approach  the 
saturation  rates. 

Experimantal  Procedure 

In  order  to  obtain  the  thermal  neutron  pulse  propagation  charac- 
teristics through  the  assembly  the  cadmium  difference  technique  was 
employed  (16) ,  using  the  cadmium  shutter  located  at  the  forward  edge 
of  the  assembly.   The  epicadmium  contribution  of  the  source  was 
thereby  eliminated. 

Normally  a  completed  neutron  wave  experiment  consur.es  days  and 
even  weeks.   During  this  tine  components  tr.ay  change  enough  to  affect 
the  data  significantly.   In  order  that  one  may  have  a  standard  of 
comparison  for  reproducibility  purposes  a  continuous  mode  run  was 
made  at  all  the  axial  positions.   These  data  can  be  acquired  in 
several  hours,  a  period  in  which  variations  are  negligible.   If 
necessary,  normalization  for  minor  variations  of  the  systems  can 
then  be  made. 

The  continuous  mode  experiment  corresponds  to  the  zero  frequency 
component  of  t^e  neutron  pulse  in  that  the  amplitude  attenuation  per 
unit  length  are  the  same  (Ic)  .   If  minor  anomalies  occur  in  the  zei'o  ■ 
frequency  component  of  the  pulsed  data,  the  dats  can  be  easily 
normalized  to  the  continuous  mode  data.   In  this  manner,  snail  changes 
due  to  day-to-day  variations  of  the  source,  e.g.,  angular  distribu- 
tions, and  counting  systems,  e.g.,  shifts  in  detector  plateaus,  can 
be  eliminated.   In  this  series  of  experiments,  no  normalization  was 
required. 
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The  continuous  node  data  which  v;ere  recorded  at  each  position  for 
both  the  cadmium  shutter  up  and  dovm  runs  included: 

(a)  run  time 

(b)  reference  detector  scaler  counts 

(c)  movable  detector  scaler  counts. 

After  the  continuous  node  experiments  were  completed,  the  pulsed 
experiments  were  conducted.  At  this  point  the  multichannel  analyzer 
which  vjas  set  up  in  accordance  with  Appendix  B  v;as  used.   Once  the 
pulsed  mode  has  been  established,  one  proceeds  to  examine  the  propaga- 
tion of  the  therraal  neutron  pulses  through  the  assembly,  employing 
the  same  cadmiuiii  difference  technique  used  in  the  continuous  mode  case- 

The  pulsed  mode  data  v;hich  were  acquired  for  each  run  included: 

(a)  run  time 

(b)  number  of  triggers 

(c)  reference  detector  scaler  counts 

(d)  movable  detector  scaler  counts 

(e)  multichannel  analyzer  printed  output 

(f)  multichannel  analyzer  binary  tape  output. 
Verification  checks  which  wera  performed  as  soon  as  the  data  were  re- 
corded are  describ'id  in  Appendix  C. 

During  each  run,  two  signals  were  continuously  monitored  with 
oscilloscopes.   One  was  the  target  current  of  the  neutron  gun.   The 
other  was  the  input  to  the  amplifiers  of  the  movable  detector  systems. 
Thus,  two  visual  presentations  vjere  available  to  assure  proper  pulsed 
operation.   A  secondary  purpose  of  monitoring  the  amplifier  input 
signal  was  to  measure  the  effective  counting  time  of  the  movable 
detector.   The  effective  counting  time  is  discussed  in  Appendix  C. 
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In  both  cadmium  covered  and  bare  pulsed  runs,  a  minumum  of 
64,000  counts  were  accumulated  in  the  peak  channel  of  the  multichannel 
analyzer.   This  sets  the  time  required  for  each  run.   Data  were  ac- 
quired every  5  cms  from  50  through  90  cms.   These  limits  are  based  on 
the  continuous  mode  analysis  which  showed  definite  deviations  from  a 
simple  exponential  decay  with  distance  for  positions  outside  this 
region.   Deviation  at  the  lower  end  was  caused  by  the  source  v;hile 
end  effects  were  significant  at  the  higher  end.   The  procedure  em- 
ployed in  determining  the  appropriate  heterogeneous  parameters  is 
outlined  in  Figure  6.   The  data  analysis  is  described  in  Chapter  IV. 
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CHAPTER  IV 

DATA  ANALYSIS  AND  RESULTS 
General 

Data  analysis  was  accomplished  with  a  modified  version  of  the 
Moore  Moments  Code  (6)  and  a  series  of  supplemental  programs 
written  for  the  IBM-I8OO.   By  using  the  I8OO  the  direct  interface 
capability  between  the  paper  tape  data  and  the  computer  could  be 
utilized.   A  flow  chart  of  the  computer  programs  is  shown  in 
Figure  7- >  while  more  detailed  descriptions  are  given  in  Appendix  D. 

Purity  of  the  heavy  water  was  estimated  to  be  99.5  ±  -S^  for 
the  moderator  only  experiment  and  99.0  +  .2%  for  the  one  fuel  rod 
experiment.   The  99.37o  purity  was  based  on  the  theoretical  model 
and  the  experimental  data.   This  purity  yielded  the  experimentally 
observed  values  of  a  and  |  at  near  zero  frequencies,  where  better 
agreement  exists  between  theory  and  experiment,  when  a  and  t  were 
computed  from  equation  [2.12].   The  99.0'>'»  purity  is  based  on  the 
difference  of  0.5%  purity  in  sanples  from  the  two  experiments  'Jhich 
were  analyzed  by  the  Departrr.ent  of  Chemistry.   Purities  were  re- 
ported to  be  99.5  +  .05%  and  98.8  +  .05%  for  the  two  experiments  (18) 
In  both  experiments  the  heavy  water  temperature  was  21  +  2°C. 

Experim.ental  errors  in  measuring  Q!  and  |  vere  computed  based  on 
the  mean  square  deviation  of  the  data  points  from  the  computed  least 
squares  fit.   Errors  ranged  from  0.5%  to.  2.0%  and  the  latter  vss 
taken  as  the  experimental  error  for  all  determinations  of  CC   and  ^. 

Resolution  time  measurements  were  conducted  as  described  in 
Appendix  A.   The  subsequent  analysis  showed  this  method  to  be  in 
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error  and  a  discussion  of  the  error  is  also  included  in  Appendix  A. 
Dufiois  (7)  measured  the  resolution  times  of  the  two  detector  systems 
just  prior  to  this  series  of  experiments.   Resolutions  times  were 
determined  to  be  1.1|.2  and  1.77  microseconds  for  the  movable  and 
reference  detector  systems  respectively.   Since  the  components  of  the 
systems  had  not  been  changed,  these  resolutions  times  were  used  in 
this  analysis. 

As  was  mentioned  in  Chapter  III,  the  17  incl  graphite  stack 
that  had  been  located  between  the  thermalizing  and  subcritical 
assemblies  was  removed  to  obtain  higher  frequency  content  in  the  data. 
There  being  less  moderation,  the  asymptotic  region  of  the  assembly  was 
smaller.   Dunlap  (5)  found  the  region  to  begin  at  25  cms.   Without 
the  graphite  stack,  the  asymptotic  region  began  at  50  cms. 

The  asymptotic  region  was  established  by  scanning   the  data. 
The  undesirable  characteristics  of  higher  spatial  and  energy  modes  at 
axial  positions  near  the  source  and  end  effects  at  positions  far  from 
the  source  tend  to  increase  the  computed  values  of  a  and  |.   In  regions 
where  neither  effect  is  significant  minimum  values  of  Zi   and  =  are 
computed.   These  minirnum  values  vere  determined  by  computing  a  and  | 
using  a  four  poinc  scan  of  the  data.   Then  a  five  point  scan  was  used 
and  the  minima  vere  compared  with  the  previous  ones.   If  they  differed 
by  less  than  l/i  the  additional  data-  point  was  considered  to  be  in  the 
asymptotic  region.   Progressively  higher  point  scans  were  used  until 
tVie  minimal  differences  exceeded  17,.   The  daua  point  producirg  this 


1.   See  the  Special  Controls  Section  of  The  ALPHA  program  des- 
cription in  Appendix  D. 
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deviation  was  excluded.   The  YL   criterion  is  based  on  the  results  of 
Booth,  Perez  and  Hartley  (16)  who  found  a  1%  tolerance  in  the  measure- 
ment of  a   and  |.   A  consistent  asyrmtotic  range,  for  al]  frequencies 
in  both  experiroants,  v;as  found  to  be  50-75  cms.   The  end  point  was 
the  same  as  that  determined  by  Dunlap  (6) . 

In  some  cases  deviating  trends  in  the  amplitudes  and  phases  were 
obviously  present.   Deviations  at  a  given  axial  position  increased  or 
decreased  with  increased  frequency.   In  these  cases  the  data  points 
were  eliminated,  beginning  at  certain  frequencies.   The  criteria  for 
elimination  were: 

a)  the  obvious  deviating  trend  at  maximum  frequency 

2 

b)  the  existence  of  a  0.257c.  difference  in  a  and  ^  between 

the  least  squares  computations  which  included  and 

excluded  the  data  point  in  question. 
Once  a  0.257,  deviation  was  observed  the  data  for  higher  frequencies  at 
the  given  position  were  also  eliminated. 

In  both  experiments  an  expansion  of  the  real  and  imaginary  com- 

2 
ponents  of  the  experimental  p   in  powers  of  ^(rad/sec)  was  made.   The 


expansions  are 


2    2  9      4 


2.      Deviations   greater    than   0.25%  produced   significant   dis- 
continuities   in   :i  and    j   when   suhseq-jent   data    points   were    eliminated. 


40 


and 

2a|  =  P^Oi     +   P  m^. 

The  P's  were  determined  by  using  the  UFNLLS  Code  (20).   Data  above 
500  cps  were  excluded  in  the  polynomial  fitting  because  statistical 
fluctuations  were  deemed  too  significant. 

Perez,  Ghanian  and  Dunlap  (19)  have  found  that  greater  sensitivity 

of  the  dispersion  law  to  thermalization  and  diffusion  parameters 

2 
exists  when  it  is  considered  in  the  p  plane  rather  than  in  the  p  plane, 

No  analysis  will  be  conducted  in  this  manner  because  the  experimental 

2 
errors  appear  to  be  more  significanc  when  the  p   dispersion  law  is 

2 
considered.   The  dispersion  laws  in  the  p  plane  will  be  shown  to 

aid  any  future  studies  that  may  be  conducted  along  these  lines. 

The  remainder  of  this  chapter  will  be  divided  into  four  parts. 

The  first  part  will  deal  with  the  continuous  mode  analysis.   The  pure 

moderator  and  one  fuel  rod  pulsed  experiments  will  then  be  analyzed 

and  finally  a  lattice  of  fuel  rods  will  be  considered. 

Continuous  Mode  Analysis 

The  continuous  mode  data  serve   three  purposes: 

1)  it  offers  a  standard  for  comparison  with  the  zero 

frequency  component  of  the  thermal  neutron  pulse 

2)  it  provides  a  standard  for  normalization  of  minor 

anomalies  in  the  pulsed  data  due  to  source  or  count- 
ing system  variations 
5)   it  establishes  a  standard  for  cor^.parlson  with  previous 
and  subsequent  experiments  . 
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It  becomes  apparent  that  the  continuous  mode  data  are  the  key  to  the 
reproducibility  of  the  experiment. 

At  each  data  position  the  thermal  ratio,  TR,  is  computed  from 
the  relation 


TR  =  I   mov   I      -  I  mov 
ref   J  bare  \  ref 


cd 


where 

N  =   the   resolution  corrected   scaler   counts 

mov  =   the   movable    detector   system 

ref  =   the   reference   detector   system 

bare  =    the   bare   run   data 

cd  =   the   cadmium  run   data. 

TR,  as  a  function  of  distance  was  fitted  to  an  exponential  function 
using  the  least  squares  technique.  The  logarithr.uc  decay  constants 
for   each  experiment   are    given    in   Table   5- 

TABLE    3. 
CONTINUOUS   MODE    DECAY  CONSTANTS 

Experiment  Asymptotic   Range  Alpha    (cm     ) 

Heavy  Water  50-75   eras  -   0.0571 

Heavy  Water-  50-75   cms  -   0.0385 

One  Fuel   Rod 

As   will   be    seen,    there   was   no   normalization    required    for   anomalies 
in   the    data.      In   this   case    the   continuous   mode   data   are   used  only   for 
comparative   purposes. 


k2 


Heavy  Vfater  Pulsed  Experiment 

The  thermal  neutron  pulses  obtained  in  the  heavy  water  only 
experiment  are  shown  in  Figure  8.   The  last  512  channels  contain  few 
or  no  counts.   These  are  insignificant  compared  to  the  peak  counts 
of  approximately  lv5>000  counts  and  are  not  plotted.   The  pulses  are 
normalized  to  the  50  era  pulse.   Each  channel  represents  50  micro- 
seconds.  These  pulses  are  the  ones  that  were  numerically  Fourier 
Transformed. 

The  Fourier  Transformed  fluxes  and  the  least  squares  fit  to  the 
amplitudes  and  phases  are  shown  in  Figures  9.  and  10.   The  maximum 
deviations  of  points  included  in  the  least  squares  fit  from  the  least 
squares  fit  are  k7,   in  a  and  17,  in  e.   Theoretical  values  of  a   and  | 
computed  from  equation  [2.12]  are  listed  in  Table  !+.  while  the  result- 
ing experimental  values  are  listed  in  Table  5.   Theoretical  and  ex- 
perir.ental  values  of  a  are  illustrated  In  Figure  11.   A  maximum  devia- 
tion of  2.1%  occurs  at  300  cps.   The  smaller  deviations  at  higher 
frequencies  is  fortuitous.   Statistical  deviations  at  these  higher 
frequencies  are  too  significant  to  allow  valid  analyses. 

An  illustration  of  :,  similar  to  the  one  of  D:,  is  given  in 
Figure  12.   In  this  case  Che  experi-eutal  data  agree  with  theory 
within  1.5%. 

The  heavy  water  dispersion  laws  are  shown  in  Figures  13.  and  \h. 
The  former  shcv;s  tr,e  experimental  and  theoretical  dispersion  laws  of 
this  work  and  the  latter  shows  a  co.^parison  of  this  work  with  that 
of  Dunlap  (c) .   Both  conparisoas  are  reasonably  good.   The  better  high 
frequency  content  of  these  data  over  that  of  Dunlap' s  (c)  is  demon- 
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THEORETICAL 
FOR 


VALUES    OF    ALPHA 
HEAVY    WATER 


AMD    Xi 


FREQUEIJCY(CPS) 

ALPHA(1/Ci1) 

0. 

-0.0375 

10. 

-0.0377 

20. 

-0.033U 

30. 

-0.0593 

uo. 

-0.0ii05 

50. 

-0.01*19 

60. 

-0.0t433 

70. 

-0.0!*ti7 

80. 

-0.01*62 

90. 

-0. 01*77 

100. 

-0.01*91 

110. 

-0.0505 

120. 

-0.0520 

130. 

-0.0531* 

lUO. 

-0.051*7 

150. 

-0.0561 

150. 

-0.0571* 

170. 

-0.0537 

180. 

-O.OGOO 

190. 

-0.0512 

200. 

-0.0521* 

210. 

-0.0536 

220. 

-0.05^3 

250. 

-0.0550 

2li0. 

-0.0571 

250. 

-0.0G23 

250. 

-0.059^ 

270. 

-0.0705 

2G0. 

-0.0715 

290. 

-0.0725 

3  00. 

-0.0737 

310. 

-  0  .  0  7 1*  7 

320. 

-0.075S 

330. 

-0.0  75  3 

3U0. 

-0.0773 

350. 

-0.C7SS 

XI  (1/CM) 

0.0000 
001*1 


-0, 
-0, 
-0, 
-0, 
-0. 

-0 
-0 
-0 
-0 
-0 


■0 
•0 
•0 
•0 


■0, 
•0, 
■0, 

■0, 
■0 


0081 
0118 
0153 
0186 

0215 
021*3 
0259 
0291* 
0317 


•0.0339 
■0.0359 
■0.0379 
■0.0598 
■0.01*15 


01*31+ 
01*51 
01*57 
01*83 


-0.01*99 


0511* 
0520 
0  51*3 
0557 
0570 


-0.05^1* 
-0.0597 
-0.0509 
-0.0522 

-0.0531* 


051*6 
0653 
0570 
0581 


■0.0592 
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TA3LE    5. 

EXPERIflEIlTAL    VALUES    OF    ALPHA    AfJD    X 
FOR    HEAVY    WATER 


FREQUEriCY(CPS) 

ALPHA  (1/Ct 

0. 

-0.037G 

10. 

-0.0378 

20. 

-0.0383 

30. 

-0.0391 

UO. 

-0.0U02 

50. 

-o.omu 

60. 

-0.0423 

70. 

-0.01+43 

80. 

-0.0458 

90. 

-0.0U73 

100. 

-0.0483 

110. 

-0.0503 

120. 

-0.0519 

130. 

-0.0534 

lUO. 

-0.0549 

150. 

-0.0564 

160. 

-0.0573 

170. 

-0.0593 

180. 

-0.0506 

190. 

-0.0619 

200. 

-0.0631 

210. 

-0.0543 

220. 

-0.0555 

230. 

-0.0569 

2tt0. 

-0.0682 

250. 

-0.0695 

260. 

-0.0703 

270. 

-0.0713 

230. 

-0.0730 

290. 

-0.0741 

300. 

-0.0753 

310. 

-0.0762 

320. 

-0.0759 

330. 

-0.0775 

340. 

-0.0782 

350. 

-0.0785 

XI  (1/CM) 

0.0000 
-0.0040 
-0.0079 
-0.0117 
-0.0153 
-0.0187 

-0.0218 
-0.0247 
-0.0275 
-0.0300 
-0.0324 

-0.0347 
-0.0358 
-0.0388 
-0.0407 
-0.0425 

-0.0442 
-0.0458 
-0.0473 
-0.0483 
-0.0503 

-0.051S 
-0.0533 
-0.0547 
-0.0551 
-0.0572 

-0.0532 
-0.0592 
-0.0503 
-0. 0513 
-0.0623 

-0.0631 
-0.0541 
-0.9552 
-0.0654 
-0.0677 
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strated  by  the  mere  consistent  results  in  the  2^10-510  cps  range. 
Results  of  the  two  experiments  differ  by  3.2%  in  a,  a  phenomenon 

attributable  to  a  difference  in  heavy  water  purity. 

2 
The  imaginary  and  real  components  of  p  are  shown  in  Figures  15. 

and  16.  respectively.   The  polynomial  fits  are  included  in  each  plot. 

Values  obtained  for  the  expansion  coefficients  are  listed  in  Table  6. 

Resultant  values  for  the  therraalization  and  diffusion  parameters  are 

included.   Values  obtained  by  Dunlap  (6)  are  also  listed.   P  and  P 

agree  within  7.1%  and  1.2%  respectively  but  the  other  coefficients 

differ  radically  in  magnitude  and  sign.   An  examination  of  Figure 

2 
16.  shows  that  the  real  component  of  p   is  erratic,  exhibiting 

oscillations  and  increasing  less  rapidly  than  Dunlap's  results  (6). 

2 
The  experimental  values  of  the  real  and  imaginary  components  of  p  ' 

are  listed  in  Table  7-   The  experimental  and  theoretical  dispersion 

2 
laws  for  heavy  water  in  the  p  plane  are  shown  in  Figure  I7.   The 

vertical  line  result  of  Age-Diffusion  Theory  is  due  to  the  neglect 
of  thermalization  effects.   The  experimental  dispersion  law  obtained 
by  Dunlap  (6)  is  also  shown.   The  theoretical  dispersion  law  obtained 
by  Perez,  Ohanian  and  Dunlap  (19)  is  also  included.   The  latter  in- 
cludes thermalization  effects.   Apart  from  the  shift  in  1   due  to  a 
difference  in  heavy  water  purity,  the  two  sets  agree  relatively  wall 
at  low  frequencies  but  pronounced  differences  occur  at  high  fre- 
quencies.  The  anomalous  behavior  of  the  heavy  water  data  in  this 
experiment  is  more  apparent  here. 

In  the  heterogeneous  analysis  the  dispersion  law  in  the  p  plane 
is  required.   This  agrees  reasonably  well  with  both  theory  and 
Dunlap's  data  (c) . 
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TABLE  6. 

EXPANSION  COEFFICIENTS  OF  RHO  SQUARED  AND  THERiVAL  I ZAT I  ON  AND  DIFFUSION 

PARAMETERS  FOR  HEAVY  V-'ATER 


QUANTITY  EXPERI:;ENTAL       ■•     EXPERIMENTAL  VALUE 

VALUE  OF  DUNlAP  (16,  19) 

pQ  [cm"'^]  1.3746  X  I O"^  1.4713  x  I O'^ 

2 
3 


P,  [cm~^  secj  5.0613  x  I0~^  5.0063  x  I0~^ 

P^  Ccm"^  sec^]         -3.1374  x  I 0~  "         3.9684  x  I O"  " 

P,  [cm""  sec^]  8.443  x  I O" '  ^  1.8025  x  lO"''' 


P4  Ccn"^  sec^J  4.3259  x  I0~'^  2.5045  x  lO"'^ 
a^^  [sec"' J                20  20 

D  [cm^  sec"' J  1.976  x  10^  1.996  x  10^ 

Cq  [cm'*  sec"' J  ***  3.73  x  10^ 


*** 


Lack  of  high  frequency  content  of  the  heavy  water  data  prevented 

a  suitable  deter-ni  nation  of  C,. . 

0 
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TABLE    7. 

EXPERIIlEfJTAL    VALUES    OF    REAL    AND 

IflAGIMARY    COMPOfiEf.'TS    OF    RflO    SQUARED 

FOR  HEAVY  WATER 


OMEGA(RAD/Sr:C) 

OMEGA  SQUARED 

REAL 

IMAGINARY 

0.000 

0.0 

0.0014139 

0.0000000 

62.831 

3947.8 
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One  Fuel  Rod  Pulred  Experiment 

In  Figure  19.  the  thermal  neutron  pulses  are  shown  lor  the  various 
axial  positions.  Again,  the  last  512  channels  are  not  shown  and  each 
channel  represents  50  microseconds.  All  pulses  are  normalized  to  the 
pulse  at  50  cm. 

Figures  19- and  20  show  the  Fourier  transformed  fluxes.   At  each 
position  considered  in  the  final  analysis,  the  amplitude  and  phase  of 
the  frequency  components  are  shown.   Included  are  the  least  squares 
fits  to  the  data.   Deviations  between  the  least  squares  fit  and  the 
data  points  considered  in  the  least  square  fit  are  scall,  the  largest 
being  less  than  57,  in  amplitude  and  less  than  17,  in  phase.   The 
theoretical  results  for  Zi.   and  |  are  given  in  Table  8.  and  the  experi- 
mental results  are  given  in  Table  9.   The  results  of  both  are  illus- 
trated in  Figures  21. and  22.   The  deviations  at  high  frequencies 
caused  by  the  neglect  of  thertr^liza  tion  effects,  are  seen  in  the 
illustrations.   T^.eery  under  estlT^tes  Ci   and  over  estimates  ^  at 
high  frequencies.   Maximum  deviations,  however,  are  less  than  5.5"  in 
OL   and  6.3%  in  |.   These  deviations  occur  at  310  cps  and  520  cps 
respectively.   The  closer  agreement  between  theory  and  experiment  at 
higher  frequencies  is  again  fortuitous.   In  the  subsequent  analyses 
data  above  500  cps,  where  maximum  deviations  were  5.1%  in  a  and 
l4..9%  i::  5,  were  excluded. 

Both  the  experirr.eatal  and  theoretical  dispersion  laws  are  sho-.'n 
in  Figure  23.   The  theoretical  dispersion  law  is  obtained  from 
equation  [2.12]  with  y  =   0.25.   This  value  of  y   gave  agreement  at 
zero  frequency  with  the  experimental  data  and  the  resulting 
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TA3LE    8. 

THEORETICAL    VALUES    UF    ALPHA    AMD    / I 
FOR    HtAVY    V/ATER-0:4E    FUEL    ROD    SYSTEM 
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FREQUEMCY(CPS) 

ALPHA(1/Cf1) 

0. 

-0.0379 

10. 

-0.0381 

20. 

-0.C3S7 

30. 

-0.0307 

1*0. 

-0. ouos 

50. 

-0.0421 

60. 

-O.Oi+3'j 

70. 

-0.04  49 
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250. 
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260. 
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270. 
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300. 

-0.0737 

310. 
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-0.0757 

330. 
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31*0. 

-0.077S 

350. 

-0.0783 

XI(1/CM) 
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-0.0078 

-0.0115 

-0.0149 

-0.0181 

-0.0211 

-0.0238 

-0.0264 

-0.0288 

-0.0311 

-0.0332 

-0.0353 

-0.0373 

-0.0391 

-0.0403 

-0.0427 

-0.0444 

-0.04G0 

-0.0475 

-0.0491 

-0.0505 

-0.0520 

-0.0534 

-0.0543 

-0.C552 

-0.0575 

-0.0533 

-0. 0600 

-0.0513 

-0.0525 

-0.0637 

-0.0648 

-0. 0650 

■0.0G71 

•0.0oS2 

TABLE    9 . 

EXPERIMENTAL    VALUES    OF    ALPHA    AflD    XI 
FOR    HEAVY    WATER-ONE    FUEL    ROD    SYSTEM 
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FREQUENCY(CPG) 

ALPHA  (l/Cf,) 
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dispersion  law  ia  shovm.   Again,  the  deviation  at  high  frequency  is 
seen.   Since  the  determination  of  7  is  based  on  the  zero  frequency 
data  where  theory  and  experiment  better  agree,  the  high  frequency 
deviation  has  no  significant  effect  on  the  experimental  value. 

The  sensitivity  of  the  fit  to  7  was  determined  by  computing  the 
dispersion  lav;  for  several  values  of  7.   It  was  found  that 

6  7  =  22  5  a' 
where  a'    is  the  zsro  frequency  a.   Assuming  that  a  is  accurate  to 
+  0.0008  (2.1%  error)  this  sensitivity  allows  the  determination  of  7 
to  +  0.02.   The  feasibility  of  using  the  neutron  wave  technique  to 

measure  7  based  on  the  dispersion  lav;  is  demonstrated. 

2 
Attention  is  now  turned  to  the  polynomial  expansion  of  p   in 

powers  o).   The  expansion  coefficients  obt.-ained  from  the  UFNLLS  Code  (20) 

are  listed  in  Table  IQ. 


table;  10. 

EXPANSION   COEFFICIENTS   C7    RKO   SQUARED 
FOR  HEAVY  VATER   -   ONE   FUEL   ROD   SYSTEM 

COEFFICIENTS  LTNTLS    RESULT 

P        Lcm'^l  l.i+28$  X    lo"^ 

0     ■■ 


P    [cm  "secj 


k.Sy^D   X  10 


P^  [cm'^sec^]  6.6175  X  lo"'''^ 

P,  [ca'^sec^l  1.2051  X  lo"'"'''^ 

P,  [cm'^ssc^]  5.3555  X  lo"''"'' 

■4' 
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Figures  2k.  and  25.  sliow  the  experimental  data  and  the  least  squares  fit 

2 

for  the  imaginary  and  real  parts  of  p  respectively.   The  2a|  fit  has 

a  maximum  deviation  of  1.57,  from  the  experimental  data  at  500  cps  while 

2    2 

O     -    I      has  a  maximum  deviation  of  2.57o  at  500  cps.   Experimental 

2    2 
results  for  2a;  and  a     -    ^     are  listed  in  Table  11. 

In  Chapter  II  P  and  P  were  related  to  7.   Using  these  relations, 

equations  [2.13]  and  [2.14],  7  was  computed  to  be  0.29  in  terms  of  P 

and  0.51  in  terms  of  P^ .   The  difference  of  these  values  from  the  0.25 

value  obtained  in  the  dispersion  law  fit  are  now  discussed  individually. 

The  P  value  of  7  is  based  principally  on  zero  frequency  data,  ? 

2    2 
being  the  zero  frequency  intercept  of  j(  -  |  ,   It  follows  that  fairly 

good  agreement  should  exist  between  the  two  values  of  7,  0.25  and  0.29, 
both  being  determined  iron  zero  frequency  data.   Nonetheless,  a  l6% 
differenc.  exists  in  7.   This  discrepancy  is  explained  by  recalling 
that  P  ,  the  zero  frequency  intercept,  was  a  computed  parameter  in  the 
least  squares  fit  Co  the  e.xperitr.ental  data.   This  value  is  less  than 
the  observed  intercept,  a  phenomenon  which  is  caused  by   the  minimum 
occurring  in  the  r^  -  l'^  data  at  approximately  50  cps.   The  minimum 
has  yet  to  be  explained,  there  beir.g  no  established  theoretical  basis 
for  its  exister^ce .   It  may  be  che  result  of  scrr.e  sys  teniatic  error.   Con- 
sequentlv,  the  experimental  value  of  ?   is  in  question.   It  is  noted 
that  if  the  value  of  the  observed  intercept  is  used,  7  is  com.puted  to 
be  0.214.  and  this  agrees  within  k%   with  the  dispersion  law  result. 

The  greater  deviation  of  the  ?  value  of  7  has  two  causes.   First, 
the  relation  derived  for  7  in  terms  of  ?   is  dependent  on  the  full  fre- 
quency range,  ?  being  the  first  derivative  of  the  2:^^  function.   A 
more  accurate  model  would  predict  a  lower  value  of  P  .   A  lower  value 
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OMEGA(RAU/SEC)  OdEGA  SQUARED 


0 
62 
125 
188 
251.327 
314.159 


000 
831 
G53 
495 


0, 
3947, 
15791, 
35530, 
63165, 
98695.8 


REAL 

00144GG 

0014455 

0  014  415 

00143G3 

0014315 

0014232 


MAG  I  NARY 
.0000000 
.0002949 
.0005913 
.000S90C 
.0011934 
.0014990 


375.990 
439. G22 
5  0  2.554 
5  6  5.436 
623.318 


142122 
193443, 
252651 
319774 


394783.4 
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0 
0 

0 
0, 


0014231 
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0 

0, 
0, 
0, 
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0018065 
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0  0  2  4  2  41 
0027320 
0030399 


691.149 
753.931 
315.313 
3  79.54  5 
942.476 


477587, 

563488, 
557133, 
7  7  3  7  7  5, 
838262, 


0 
0 
0 
0 
0 


0014643 
0014774 
0014921 
0015112 
0015543 


0 
0 
0 

0 
0, 


0033485 
0036562 
0039552 
0042771 
0045890 


1005 
1053 
1130 
1193 


303 
140 
972 

8  04 


1256.636 


1010645 
1140924 
1279093 
14  2  515  3 
1579133 


0, 

0, 

0 

0, 

0, 


0015676 
0015016 
0015230 
0015544 
0016773 


0 
0 
0 
0 


0043930 
0051957 
0054947 
0057953 


0.0061051 


1319 
1332 

1445 
1507 


458 

2  99 
131 
96  3 


1570.795 


1749994. 7 
1910751.7 
2  0  3  3  4  0  4.1 
2273952.7 
2467396.0 


0017161 
0017504 
0013203 
0  013  9  2  7 
0019571 


0, 
0, 
Q, 
0, 
0, 


0  0  5  4  2  3  4 
0067373 
0070574 
00  735  53 
0076470 


1533 
1696, 
1759 
1322 


62  6 
45C 
2  90 
122 


1384.953 


266: 
2377971.1 
30  9  5101.9 
3  3  2  012  3.4 
3553050.  7 


0 
0 

0 

n 

0 


0020113 
00203  63 
0  0  215  91 
0022690 
0024522 


0 
0, 

0, 

0 

0, 


0079417 
0CS2553 
0035450 
0083723 
0091399 


1947, 
2010 
2073 
213  5.231 
2199.113 


735 
517 
44  9 


3793853.0 
4  0  4  2  5  31.5 
4299190.4 

45635  9  6.8 
4335095.9 


0.0025514 
0.0026108 
O.C025451 
0.0024430 
0  .  C  0  2  1 5  2  9 


0.009  4  011 
0.0095803 
0.0093191 
0,0100510 
0.0104173 


72 

of  P^  produces  a  lower  value  of  7.   7  computed  from  the  P  relation 
was  found  to  be  very  sensitive  to  the  value  used  for  v   .  2,kQ0   m/sec 
and  2,1|90  m/sec  yielded  results  of  7  =  0.2?  and  7  =  0,22  respectively, 
The  2,470  m/sec  is  the  average  therm.al  neutron  velocity  in  heavy  water, 
Insertion  of  the  fuel  rod  would  tend  to  increase  the  average  velocity 
because  of  preferential  absorption  of  lower  energy  neutrons  and  pro- 
duction of  fast  neutrons  in  the  fuel  rod.   Consequently,  the  correct 
value  of  v^^   should  be  somewhat  larger  than  that  for  the  moderator, 
resulting  in  a  more  consistent  value  of  7. 

For  the  present,  7  =  0.25  +  .02  obtained  from  the  dispersion  law, 
is  considered  the  most  accurate  value.   Assuming  that  Diffusion  Theory 
holds  in  the  rod,  7  was  computed  to  be  0.3227  by  Cain  (15). 

The  experimental  and  theoretical  dispersion  laws  in  the  p^  plane 
are  shown  in  Figure  26.   The  high  frequency  difference  is  amplified. 
The  next  analysis  deals  with  the  intersection  of  the  dispersion 
laws  of  the  heavy  water  ouly  and  heavy  water-one  fuel  rod  systems. 
This  analysis  is  hampered  somewhat  by  the  difference  in  heavy  water 
purity  that  existed  between  the  two  sets  of  data.   This  difference 
was  taken  into  account  by  reducing  the  cc  of  the  one  fuel  rod  data 
0.00157  to  compensate  for  absorptions  of  the  additional  0.5%  light 
water.   This  correction  was  obtained  by  computing  the  zero  frequency 
a  at  the  two  concentrations  of  heavy  water  from  equation  [2.12].   The 
intersection  is  then  based  on  a  99.5%  heavy  water  purity.   The  dis- 
persion law  intersection  is  shown  in  Figure  27.  where  both  the  experi- 
mental and  theoretical  intersections  are  shown.   Tneore tical ly  the 
intersection  occurs  at  310. Ii  cps,  a  value  nrjch  hieher  than  the 
experimental  intersection  at  158  cos  which  was  obtained  by  plotting 
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the  observed  dispersion  laws  for  the  heavy  water  only  and  heavy 
water-one  fuel  rod  systems.   This  discrepancy  is  due  to  the  experi- 
mental error  and  the  inadequacy  of  the  Age-Diffusion  Model.   The 
dominant  experimental  error  is  caused  by  the  uncertainty  of  the  heavy 
water  purity  for  it  was  found  that 

&o  -x,  7,000  5C 
c    ' 

where  C  is  the  correction  term  applied  for  the  purity  difference  be- 
tween the  two  experiments.   C  is  approximate  and  introduces  signifi- 
cant deviations.   More  accurate  results  would  have  been  obtained  if 
the  experiments  had  been  conducted  at  the  sane  heavy  water  purity. 
On  the  other  hand,  an  improved  theoretical  model  would  predict  a 
smaller  value  of  03   thereby  producing  better  agreement  between  the 

theoretical  and  experimental  intersection.   The  intersection  ia  the 

2 
p  plane  is  shown  in  Figure  23.   Uncertainty  in  the  heavy  water 

results  obscure  this  analysis  but  a  more  distinct  intersection  may 

occur. 

The  combination  of  theoretical  and  experimental  uncertainties 

prevent  the  accurate  determination  of  heterogeneous  parameters  based 

on  the  intersection  of  the  dispersion  laws.   This  is  readily  apparent 

in  Table  12.  where  the  experiniental  and  predicted  parameters  are  given. 

A  more  accurate  experiment  and  theoretical  model  are  needed  for 

proper  analysis  of  the  dispersion  law  intersection. 

A  Lattice  of  Fuel  Rods 

7  has  been  measured  for  a  one  fuel  rod  systeni.  Will  the  value  of 

7  change  when  a  lattice  of  fuel  rods  is  considered?  It  is  not  possible 

at  this  time  to  answer  this  question  resolutely  but  che  work  of 
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TABLE    12. 

HETEROGENEOUS    PARAMETERS   BASED  ON 
INTERSECTION  OF   DISPERSION  LAWS 

VARI.-vBLE  EXPERIMENTAL  VALUE  _  PREDICTED  VALUE 

a  [cm      ] 
lJcm-1] 

03    [cps] 

CD    [cps] 

r      ^^ 
P 


-0.03515 

-0.07i+78^ 

-0.03855 

-0.06575^ 

150.8 

510.0^ 

120.9 

510.4^ 

552.5 

118.8^ 

>.  1 

2 
.999 

1  These   values   vere   obtained   from   the  Age-Diffusion  Model. 

2  These   values   v;ere   obtained   froa  N.    J.    Diaz   by  private   corxiuni- 

cation. 
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Dunlap  (6)  does  offer  a   clue  to  the  answar.   His  study  included  the 

calculation  of  cell  parameters  of  a  lattice  of  fuel  rods  of  the 
same  type  as  the  one  considered  in  this  experiment.   The  calculated 
parameters  yielded  a  theoretical  dispersion  law  which  agreed  reasonably 
well  v.'ith  the  experimental  one.   The  agreement  leads  one  to  believe 
that  the  calculated  and  true  values  of  the  parameters  are  nearly  the 
same.   The  value  of  L  for  the  lattice  will  now  be  considered. 

L  of  the  lattice  is  related  to  y   and  L^  of  the  moderator  (15) 
and  is  very  sensitive  to  y    in  that 

6L^  ^  1.6  X  10^  57. 

Using  Dunlap's  data  and  7  =  0.25,  L^  was  computed  to  be  92.2  cm^. 
This  compares  within  0.6%  of  the  91.63  cm^  value  computed  by  Dunlap 
(6).   The  close  agreement  in  L^  and  the  agreement  obtained  by 
Dunlap  between  the  theoretical  and  experimental  dispersion  law  in- 
dicate that  the  7  measured  for  one  fuel  rod  may  be  applicable  to 
a  simple  lattice  of  fuel  rods. 


CH/JTER  V 
CONCLUSIONS  AIJD  RECOMivffiNDATIONS  FOR  FUTURE  WORIC 

Conclusions 

The  neutron  vave  technique  has  been  used  to  study  a  heterogeneous 
system.  7  has  been  determined  using  the  experinental  and  theoretical 
dispei'sion  lavs  and  the  expansion  coefficients  of  0   in  a  povrer  series 
of  (iw).  The  intersection  of  the  experinental  dispersion  lavs  of  the 
heavy  vater  only  and  the  heavy  water-one  fuel  rod  syster-.s  nade  possible 
the  determination  of  tvo  other  heterogeneous  paramaters.   Heterogeneous 
parameters  have  been  extracted  froni  the  neutron  vave  data.   These  para- 
meters are  sufficiently  sensitive  to  neutron  vave  paraneters  to  allov 
their  determination  by  this  technique.   Accuracy  of  this  technicue  has 
been  partially  established  in  that  v  can  be  deteniiined  through  the  dis- 
persion lav  fit  to  tvo  significant  digits.   Further  investigation  of 
accuracy  vas  precluded  because  of  the  theoretical  xodel  and/or  experimental 
phenomena. 

The  extension  of  the  results  obtained  for  v  tc  a  lattice  has  ce=n 
investigated  for  the  case  cf  a  simple  lattice.   The  diffusion  length 
squared  for  the  lattice  was  computed  based  on  the  meas'ired  value  of  y 
and  it  was  found  to  agree  v&vj   well  with  previous  results.   The  fact  that 
the  correct  value  of  the  diffusion  length  squared  can  be  computed  from 
the  measured  vaZ.ue  of  y    is  verj'  significant. 

Tne  experimental  technicue  has  been  tested  under  the  most  adverse 
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conditions  in  that  a  rdnimim  of  fuel  vas  used.  V.^ith  larger  amounts 
of  natural  uranium  or  with  enriched  uranium  more  pronounced  character- 
istics vill  be  exhibited.   In  these  systems  improved  results  should  be 
obtained. 

Recommendations 
The  most  important  recammendation  is  the  development  of  a  more 
accurate  theoretical  model.   Better  accuracy  at  high  frequencies  would 
improve  the  technique,  predicting  more  accurate  values  for  the  appropri- 
ate parameters. 

Future  studies  should  be  conducted  to  determine  interaction  and 
shadowing  effects  on  the  heterogeneous  parameters.  These  studies  can 
be  accomplished  by  a  study  equivalent  to  this  one  where  two,  three  or 
more  fuel  rods  are  inserted  at  various  locations  in  the  assembly.  An 
interesting  and  useful  study  would  be  the  deteiroination  of  heterogenous 
parameters  as  a  function  of  the  number  of  fuel  rods  in  the  asser.blv. 
Still  another  approach  weald  be  the  study  of  the  heterogeneous  para- 
meters as  a  function  of  the  center  line  separation  of  the  fuel  rods. 

Tne  analysis  of  fu-.ure  experiments  may  be  modified  in  the  follow- 
ing way.   Let  the  parametric  fit  in  7  of  the  theoretical  to  the  experi- 
mental disp-r.-icn  lavs  ce  the  method  cf  determine  7.   Then,  in  the 
evaluation  of  the  expansion  cceff  ic5  er  ts  ,e.i?  .  (2.13)  and  (2.lU),  y 
is  known  and  two  other  parameters  can  be  determined. 

The  dispersion  lavs  intersection  was  not  studied  in  the  p^  plane. 
It  may  be  worthwhile  to  do  so  in  the  future.   It  appears  that  a  more 
distinct  intersection  may  occur. 

Remaining  reccrimeniations  are  concerrea  with  the  experimental 
procedure.  An  accurate  and  reliable  method  of  deterrdng  the  resolution 
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time  of  the  counting  system  is  needed.   The  method  described  in  Appen- 
dix A  may  serve  as  a  starting  point  for  the  development  of  a  suitable 
method.   It  is  important  to  determine  the  resolution  time  accurately 
to  0.01  seconds  because  the  experimental  value  of  o.  and  E,   are  dependent 
on  the  resolution  time  used.   A  more  accurate  resolution  tine  yields 
a  more  accurate  experimental  values  of  a  and  C  and  the  other  neutron 
wave  parameters. 

Determination  of  heavy  water  purity  proved  to  be  a  difficulty  in 
this  experiment.   Should  heavy  water  or  some  other  moderator  with  vari- 
able properties  be  used  in  future  experiments,  the  physical  properties 
must  be  accurately  determined  and  verified.   Efforts  should  be  made  to 
conduct  the  moderator  only  and  moderator-fuel  experiments  under  the 
same  conditions'  when  such  moderators  are  used.   The  error  in  determin- 
ing the  physical  properties  must  be  minimized  to  improve  the  validity 
of  the  data  analysis. 


APPENDIX  A 
t-ffiASURE-^IZI-iT  OF  RZ30LUTI0IJ  TIMS  USING  A  PULSED 
NEUTROII  SOURCE 
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The  method  presented  here  for  determination  of  the  resolution 
time  was  found  to  he  inadequate  when  perforraed  in  the  manner  described 
herein.   The  method,  however,  is  believed  to  be  sound  and  with  one 
improveraent  good  results  should  be  obtained. 

The  difficulty  concerns  the  rectangular  shape  of  the  neutron 
p\ilse  at  the  detector.   The  asymptotic  decay  constant  of  the  therraalizing 
apparatus  was  such  that  the  individual  pulses  at  the  detector  are  not 
sharp  and  the  decay  time  consu:?ies  a  great  percentage  of  the  tine  between 
pulses.   Consequently,  an  "effective"  time  between  pulses  exists  rather 
than  the  time  between  the  end  of  one  pulse  and  the  beginning  of  the 
subsequent  one  which  was  measured  at  the  source.  An  improved  procedure 
to  overcoiiie  this  difficulty  is  discussed  at  the  end  of  this  Appendix. 

Resolution  time  is  defined  as  the  minimum  time  between  successive 
neutron  interaction  events  such  that  both  events  are  registered  by  the 
counting  system.  The  definition  is  not  changed  if,  instead  of  a  single 
neutron,  a  burst  of  neutrons,  as  shown  in  Figure  29.  is  used.  Assuming 
that  the  neutron  bursts  from  the  neutron  generator  are: 

(a)  the  saiie 
and 

(b)  near  a  square  wave  at  the  detector 

a  method  for  deter:nining  the  resolution  uime  of  the  counting  syster.s 
was  developed  which  used  the  neutron  generator.  The  major  advantage 
of  this  method  is  that  the  resolution  time  of  trie  counting  system  in 
its  operating  condiLion  is  determined. 

The  method  is  illustrated  graphically  in  Figure  29.   It  is  seen 
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that  when  the  time  between  the  end  of  one  burst  and  the  beginning  of 
the  next  is  less  than  the  resolution  time,  -j,  of  the  system  the  effec- 
tive number  of  neutrons  counted  per  burst  is  less,  because  some  of  the 
leading  neutrons  of  the  burst  are  not  counted.   Othei-^'ise  the  effective 
number  of  counts  per  burst  is  constant. 

To  verify  assumption  (a)  the  target  current  was  monitored  with 
a  Type   5'-tT  Tecktronix  Oscilloscope  with  a  lAl  Plug-in  Pre/Jiiplif ier 
unit.   This  Pre.toplifier  unit  permitted  double  amplification  of  the 
input  signal  by  cascading  the  two  channels  of  the  PreAmplifier.  This 
monitoring  method  presented  a  trace  of  the  target  current  which  could 
be  monitored  accurately  enou-h  for  assurance  of  constant  target  current 
and  hence  a  constant  burst  of  neutrons. 

Since  moderation  will  spread  the  neutron  burst,  the  detector  was 
placed  as  near  the  source  as  practical.   In  this  case,  the  detector 
was  placed  in  the  reference  detector  position  shown  in  Figure  3. 
This  arrangement  was  the  optim'am  position  for  obtaining  nearly  sauare 
pulses  of  thermal  neutrons  at  the  detector. 

Tne  resolution  time  of  the  system  was  believed  to  be  of  the  order 
of  several  microseconds.   For  this  reason  data  were  taken  using  a  10 
microsecond  target  pulse  initially  separated  by  an  3  microsecond  time 
interval .   The  time  intervals  were  measured  by  observing  the  terminating 
and  initiating  signals  for  the  neutron  pu.lse3  on  an  oscilloscope. 
The  time  interval  of  separation  was  decreased  until  a  nonlinear  devia- 
tion was  observed. 

The  integral  number  of  counts  were  recorded  by  the  detector  system. 
If  the  time  interval- of  pulse  separation  is  greater  than  the  resolution 
time,  the  counts  shoiold  increase  linearly  with  decreasing  time  separa- 
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tion  (increasing  repetition  rate).   '.^Iien  the  time  separation  is  less 
than  the  resolution  tiz-.e  there  should  be  a  nonlinear  deviation  in  the 
scaler  counts.   Representative  data  are  given  in  Table  ]3.  and  illustrated 
in  Figure  30. 

In  order  to  overcome  the  difficulty  of  pulse  spreading  the  genera- 
tor should  be  withdrawn  from  the  thermalizing  assembly  and  a  saall 
a-Tiount,  for  exaniple  several  inches,  of  moderator  placed  between  the 
target  and  the  detector.  Although  the  sharpness  of  the  detector  pulse 
has  not  been  completely  verified,  a  later  experiment  indicates  that 
this  procedure  will  yield  suitable  results  for  the  resolution  time. 


TABLE  13. 
DATA  FOR  RESOLUTION  TIME  OETERM  I  .MAT  I  ON 
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TIME    BETWEEN 
BURSTS     (illCROSECONDS) 


2 
3 
k 
6 
8 


AVERAGE    COUNTS 
FOR    ONE    MliNUTE 

3ti0^792 

3  (i  '1 ,  6  5  5 

3i»3,651 

318,50  7 
301,070 
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Operation  in  the  pulsed  neutron  mode  requires  a  mutual  compata- 
bllity  between  the  current  pulse  width  at  target,  the  time  svjeep  of 
the  multichannel  analyzer  and  the  thermal  neutron  decay  constant  for 
the  system  being  pulsed.   Since  one  would  like  to  have  pulses  with  the 
ratio  of  high  frequency  to  lov;  frequency  contents  as  high  as  possible, 
a  minimal  current  pulse  width  at  target  is  dictated,   however,  one 
also  needs  to  obtain  good  statistics  so  the  minimal  pulse  width  v;ill 
be  longer  depending  on  the  amount  of  time  one  wishes  to  devote  to  the 
accumulation  of  data.   For  this  set  of  experiments  a  pulse  v;idth  on 
target  of  2.0  milliseconds  was  selected  as  the  practical  compromise 
between  these  two  conditions. 

The  multichannel  analyzer  titr.e  sweep  must  be  long  enough  to  des- 
cribe fully  the  decay  to  background  of  the  neutron  population.   This 
is  most  important  since  the  data  analysis  assumes  that  the  neutron 
pulse  decays  to  zero  in  the  last  channels  of  time.   The  buildup,  decay 
and  background  levels  must  be  recorded  by  the  m.ultichannel  analyzer. 

For  the  experiments  performed  here  the  channel  width  selected 
was  kO   microseconds.   This  yielded  a  multichannel  analyzer  time  sweep 
of  51.2  milliseconds.   The  cutoff  was  found  to  be  in  the  background 
region  in  all  cases. 

Now  that  the  multichannel  analyzer  time  sweep  is  established  a 
repetition  rate  for  the  neutron  gun  is  determined.   Proper  operation 
requires  that  a  trigger  signal  initiates  the  sweep  of  the  multichannel 
analyzer  and  operation  of  the  neutron  generator  simultaneously.  Two 
milliseconJs  later  (the  width  of  the  pulse  on  target)  the  neutron 
generator's  operation  is  terminated  by  the  final  pulse.   51.2  milli- 
seconds later  the  multichannel  analyzer  completes  its  sweep.   The 
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time  for  the  next  initiating  trigger  signal  can  be  any  time  after 

51.2  milliseconds.   It  should  be  close  to  51.2  milliseconds  in  order  to 

minimize  the  lost  time  between  cycles  of  operation. 


APPEITDIX  C 
EXPERIMENTAL   CHECKS 
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Introduction 

As  soon  as  the  data  are  recorded  certain  checks,  auxiliary  neasure- 
ments  and  computations  can  be  made  to  insure  that  the  data  are  valid. 
Each  one  is  described  and  the  reasons  for  it  are  given.  The  only 
equipment  needed  is  a  desk  calculator.   Since  they  are  extremely  simple 
it  would  be  poor  Judgement  not  to  perform  them. 

Effective  Countint^  Fraction 

The  effective  counting  fraction  is  the  ratio  of  the  detection 
system  counting  tine  in  a  cycle  to  the  total  cycle  time.  The  effective 
counting  time  is  needed  in  order  to  make  proper  resolution  corrections 
which  should  be  based  on  the  actual  count  rate  rather  than  the  count 
rate  averaged  ever  the  complete  cycle.  During  the  neutron  pulse  the 
count  rate  is  at  a  maximum  value.   But  the  pulse  only  lasts  for  approxi- 
mately 10^  of  the  cycle.   After  the  decay  of  the  pulse  the  count  rate 
is  essentially  zero.   Therefore,  an  averaged  count  rate  is  aoproximately 
10>  of  the  effeccivs  count  rate  and  resolution  corrections  based  on 
the  averaged  count  rate  vould  be  too  small.   The  effective  counting 
fraction  allows  one  to  determine  the  actual  count  rate. 

The  input  to  the  main  a.~plifier  furnishes  a  suitable  oscillo- 
scope display  for  determination  of  the  effective  counting  fraction. 
A  representative  display  is  shown  in  Figure  31.  The  counting  time  is 
the  width  of  the  pulse  packet.  The  cycle  tine  is  the  time  between 
the  leading  edges  of  two  succeeding  packets.   The  effe.^tive  counting 
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fraction  F  is  defined  by  the  equation 

F  =  time  vidth  of  pulse  packet, 
cycle  time 

Resolution  corrections  are  then  computed  by  the  relation 


K  =  ^ 


1  -  HI 
TF 

vhere 

N  =  resolution  corrected  counts 

n  =  observed  counts 

T  =  total  run  time 

T  =  resolution  time. 
The  resolution  correction  of  all  pulsed  data  was  computed  in  this 
manner. 

The  effective  counting  fraction  changes  with  position  because 
the  neutron  pulse  becomes  vider  as  distance  from  the  source  is  increased < 
Therefore,  the  fraction  must  be  measured  for  the  movable  detector  syste* 
at  each  a:cial  position.   The  reference  detector  is  fixed  so  its  fraction 
has  to  be  r.-.easared  only  onze.      Ilote  that  the  fractions  will  change  if 
the  neutron  generator  target  current  pulse  width  is  changed.   The 
effective  counting  fraction  for  the  reference  detector  w?.3  C.07  while 
the  movable  cstectcr's  fraction  varied  from  0.2U  at  z=50  cm  to  0.31 
at  2=75  cm. 

Peak  to  Background  Ratio 

A  pulse  neutron  experiment  should  be  conducted  with  a  high  peak 
to  background  ratio,  1,000/1  being  a  minimijm  ra-io.   The  backgro'^:d 
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vas  obtained  'oy  averaging  over  the  last  100  channels.  The  peah  to 
background  ratio  varied  from  65,000  and  1*  ,000  at  z-50  en  to  2U,000 
and  2,500  at  z=75  en  in  the  moderator  only  and  moderator-fuel  rod 
experiments  respectively. 

Tine  Per  Trigger 
The  run  time  represents  the  time  between  each  initial  trigger  of 
the  system.  This  number  should  be  greater  than  the  multichannel  analyzer 
time  sweep  and  equal  to  the  measured  time  between  trigger  pulses.   Some 
variation  occurs  due  to  drift  in  the  time  between  triggers  and  loss 
of  trigger  signals  in  the  electronic  system.   Variations  of  greater 
thaji  l.O;';  are  indicative  of  system  msilfunction.   In  this  series  of 
experiments  the  tine  per  trigger  was  51.63  *  .Uo  or  51. c3  -  0.5;'  sec. 
Tnennal  Ratio  of  Movable  Dector  Scaler  Counts  to 
Reference  Detector  Scaler  Counts 


The  resolution  corrected  counts  computed  as  described  in  the  second 
section  of  this  Appendix  are  used  in  this  step.  The  thermal  ratio  is 
computed  in  the  manner  as  the  continuous  mode  thermal  ratio  was  computed 
in  Charter  IV. 

The  thermal  ratio  computed  here  is  compared  with  the  one  obtained 
from  the  continuous  mode  data.   The  two  should  be  comparable,  and  as 
a  function  of  distance,  the  thermal  ratios  shcold  possess  the  ssme 
linear  slope  on  a  semi-log  plot.  The  more  significant  part  of  this 
check  is  the  comparision  of  the  slopes.   If,  at  a  given  position,  the 
thermal  ratios  do  not  reasonably  compare  (within  '-  25f^)  in  magnitude 
and  the  slopes  differ  'ay   more  than  a  few  percent,  the  run  should  be 
repeated.   The  slopes  differed  by  1.6':i   in  the  moderator  only  experiment 
atnd  by  1.3"^  in  the  moderator-one  fuel  rod  experiment.   The  thermal 
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ratios  differed  by  as  nuch  as  12% 

Ratio  of  Movable  Detector  System  Multichannel  Analyzer  Counts 

to  Scaler  Counts 
The  movable  detector  system  counts  are  measured  by  both  the  scaler 
and  the  multichannel  analyzer.  Kovever,  the  multichannel  analyzer 
counts  for  a  fraction  of  the  cycle  time  due  to  dead  tine  between  channels 
and  dead  time  between  the  last  channel  and  the  beginning  of  the  next 
cycle.  The  multichannel  analyzer  counting  time  (channel  width  times 
the  number  of  channels)  divided  by  the  time  between  triggers,  described 
above,  should  equal  the  ratio  of  multichannel  analyzer  integral  counts 
to  the  movable  detector  system  observed  scaler  counts.   Some  variation 
will  occur  because  of  the  differing  resolution  times  of  the  multi- 
channel analyzer  and  the  scaler.   However,  the  ratio  should  not  vary- 
more  than  a  few  tenths  of  a  percent.   In  this  series  of  experiments 
the  ratio  was  found  to  be  O.SlOo  *  .0020  or  O.8IO8  *  0.25)^- 


APPENDIX  D 
COMPUTER  PROGR.\MS 
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General 

A  set  of  computer  programs  to  be  used  in  the  analysis  of  neutron 
wave  propagation  data  has  been  v;ritten  for  the  IBM-I80O.   This  com- 
puter has  a  paper  tape  reader  which  allows  the  paper  tape  output  of 
the  multichannel  analyzer  to  be  read  in  directly.   The  I80O  programs, 
as  well  as  those  used  on  the  IBM  56O/5O,  are  now  described.   Each 
description  includes  a  brief  synopsis  of  the  program,  the  input/out- 
put and  any  special  control  features  utilized.   Additional  comments 
are  made  at  the  end  of  the  sections.   All  programs  are  written  in 
FORTR.\N  IV. 

An  extensive  number  of  plots  are  among  the  output  of  the  various 
programs.   The  size  of  the  axes  has  been  set  at  five  inches  along  the 
horizontal  and  eight  inches  along  the  vertical  directions.   This  can 
be  changed  by  adjusting  the  appropriate  scaling  factor  and  recompiling 
the  program.   Any  type  or  size  of  paper  can  be  used  as  long  as  one 
dimension  does  not  exceed  11  inches. 

TAPRD 
Synopsis 

This  1600  program  reads  the  paper  tape  output  of  the  multichannel 
analyzer  and  stores  the  content  of  each  channel  on  magnetic  disk. 
After  the  tape  is  read  the  channel  content  is  printed  out  for  com- 
parison with  the  printed  output  of  the  multichannel  analyzer.   If  the 
tape  has  been  read  correctly,  it  is  stored  on  disk.   Otherwise, 
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another  attempt  can  be  \nade    to  read  the  tape.   A  bare  data  tape  is  read 
first,  followed  by  a  cadmium  data  tape.   The  data  are  stored  on  disk 
in  groups  of  256  channels.   Each  group  is  assigned  a  file  number. 

Input 

The  following  data  are  read  from  the  keyboard  in  FORIL-\T  (^Ik)  : 
L  =  first  storage  file  number,  usually  an  odd  3  digit  number 
tux  =  number  of  runs,  each  run  having  a  bare  and  cadmium  data 

tape 
N  =  number  of  channels  per  data  tape. 

Output 

To  identify  the  storage  files  the  initial  file  number  is  printed 
The  channel  number  and  content  of  every  50th  channel  are  then  printed 
After  all  data  tape  sets  have  been  read,  the  next  available  file 
number  is  printed. 

Special  Controls 

PAUSE  llll_^--  This  pause  occurs  so  that  the  first  (bare)  data 
tape  can  he  readied  in  the  ca.p&  reader.  'v"aea  readied,  press  START 
to  begin  the  reading  process. 

PAUSE  a.i:.Ll:._--   This  pause  occurs  after  the  printout  of  each  data 
tape  and  allo^vs  one  to  verify  that  the  tape  has  read  correctly.   If 
an  error  has  occurred,  turn  SZN5E  S''ITCH  0  on  and  press  SI.-\RT  Co 
return  the  program  to  PAUSE  1111,   If  the  tape  has  read  correctly, 
press  START  to  store  the  data  on  disk. 

PAUSE  2222.--  This  pause  occurs  so  that  the  second  (cadmium) 
data  tape  can  b3  readied  in  the  tape  reader,  '-rnen  readied,  press 
START . 

PAUSE  i|I|.kij. .  -  -   This  pause  is  the  same  as  the  previous  PAUSE  likkk 
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except  that  turning  SENSE  SWITCH  0  on  and  pressing  START  returns 
the  program  to  PAUSE  2222. 

PAUSE  y^^^.--   I'his  pause  occurs  when  all  the  data  sets, 
specified  by  MAX,  have  been  read.   Pressing  START  returns  the  program 
to  the  initial  reading  of  the  keyboard  input  data. 

Coioraent 

Approximately  I5  seconds  are  required  to  read  a  1,0211  channel  tape, 

TOERM 
Synopsis 

This  1800  program  computes  the  corrected  and  normalized  thermal 
neutron  pulse  from  the  bare  and  cadmiu:Ti  data  stored  on  disk.   Each 
channel  of  the  bare  and  cadmium  data  is  corrected  for  resolution 
losses.   The  number  of  thermal  neutrons  per  channel  is  then  computed, 
background  is  subtracted  and  nornialization  is  performed.   Due  to  the 
limited  core  capacity  of  the  IcOO  the  computations  are  made  in  groups 
of  256  channels.   The  complete  normalized  pulse  is  stored  on  disk  in 
one  file  after  each  group  of  corzputa tions. 

Inout 


Input  data  is  read  from  cards.   One  initial  set  of  two  cards  is 
required  for  each  set  of  runs.   A  second  set  of  two  cards  is  required 
for  each  run  within  the  set.   The  initial  cards  contain  the  following 
data  in  FOR>L\T  (kElO.ii,  2E13. 7/6li+)  : 

TAU  =  resolution  time  of  movable  detector  system  in  m.icroseconds 
CTAU  =  resolution,  time  of  reference  detector  system  divided  by 

the  effective  counting  fraction  of  the  reference  detector 
system  (See  Appendix  D) 
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TI-ICCH  =  channel  width  of  multichannel  analyzer  in  microseconds 

TIICDT  =  deadtir.e  between  channels  in  microseconds 

STREF  =  observed  reference  detector  system  counts  for  bare  run 
to  which  the  data  are  normalized 

TIMER  =  run  tirre  for  which  STREF  counts  were  observed 

NCH/\N  =  number  of  channels  of  data 

NSKIP  =  print  control  parameter,  every  NSKIPth  channel  is 
printed  beginning  with  channel  ^1 

MRUNS  =  meximum  number  of  runs  in  the  set 

I.l    =  initial  file  nunber  of  the  bare  data  of  the  first  run  for 
which  the  thermal  pulse  is  to  be  computed 

L2    =  the  file  number  for  storing  the  first  thermal  pulse  com- 
puted in  the  set,  usually  a  1+  digit  number  ending  in  1 

IBKGD  =  the  nurVoer  of  final  channels  over  which  the  average 
background  is  computed. 
The  run  data  cards  contain  the  following  numbers  in  FO:^L\T  (tE13.7/3^^) 

REFl  =  observed  reference  detector  system  counts  for  bare  run 

REF2  =  observed  reference  detector  system  counts  for  cadmium  run 

TDiEl  =  time  for  bare  run 

TI>2;2  =  time  for  cad.T.ium  run 

TRIGl  =  number  of  triggers  for  bare  run 

TRIG2  =  number  of  triggers  for  cadmium  run 

NRUN   =  run  number 

NZPOS  =  axial  position 

N1^-\D  =  radial  position. 
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Output 

The  observed  bare  counts,  resolution  corrected  bare  counts, 
observed  cadmiun  counts  and  resolution  corrected  cadmium  counts,  are 
printed  out  for  every  NSKIPth  channel.   Then  the  normalized  and 
background  corrected  thermal  neutron  pulse  is  printed  out  for  every 
NSKIPth  channel.   The  last  portion  of  the  printout  is  the  background 
per  channel  of  the  bare,  cadmium  and  thermal  pulses.   After  all 
pulses  have  been  computed  the  next  available  file  number  is  printed. 

Special  Controls 

PAUSE  •^^^^.--   This  pause  occurs  when  the  thermal  neutron  pulses 
have  been  computed  for  MRUNS  runs.   Pressing  START  returns  the  pro- 
gram to  the  reading  of  the  initial  set  of  data  cards. 

Comments 

The  logic  of  the  file  numbering  system  can  now  be  explained. 
Assume  that  one  has  10  sets  of  l,021i  channel  data  tapes  and  that 
L  =  lOl  in  the  TAPRD  input.   Four  files  will  be  required  to  store  each 
data  tape  ik   X  2^6  =  1,021^.).   The  first  bare  data  tape  will  be 
store-:  in  files  =101,  #103,  =105,  and  #10?.   The  first  cadmium  data 
tapes  will  be  stored  in  files  =102,  #10^.,  ylOo,  and  -flo3.   Files  =101 
and  ^102  contain  the  bare  and  cadmium  data  for  channels  1  through  256. 
Succeeding  data  are  stored  in  the  same  fashion  in  higher  file  numbers. 
If  L2  =  1001  in  the  TOZ?!:'!  input  then  the  normalized  pulses  will  be 
stored  in  files  ;fl001  through  -^lOlO  inclusive.   Subsequent  sets  can  be 
identified  by  using   201  and  2001,   5OI  and   3OOI,  etc.,  for  Lland 
L2  respectively.   Thus,  the  leading  digit  can  bs  used  to  identify 
the  different  data  sets. 
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PLOT! 
Synopsis 

This  1800  program  plots  the  thermal  neutron  pulses  which  were 
computed  and  stored  by  THERM.   Scaling  and  labeling  is  based  on  the 
first  pulse  plotted  (the  LMIN'th  file)  which  is  usually  of  raaximura 
amplitude. 

Input 

The  following  data  are  read  from  the  keyboard  in  FOfl>L'vT  (I314)  : 
LMIN  =  file  number  of  first  pulse  to  be  plotted 
LMAX  =  file  number  of  final  pulse  to  be  plotted 
NCHAN  =  number  of  channels  to  be  plotted 
LBAD  =  file  number  of  pulses  to  be  omitted  in  plotting, 
a  maximum  of  10,  leave  blank  if  none. 

Output 

The  thermal  neutron  pulses  are  plotted,  superimposed  on  the  same 
axes . 

Special  Control 

_PAUS£  ;,  722. --This  pause  occurs  after  all  plots  are  made.   Pressing 
START  returns  the  program  to  the  reading  of  the  keyboard  input  data. 
Conments 

The  X  axis  is  labeled  at  five  points  including  zero.   To  obtain 
integer  numbers,  NCR^N  should  be  divisible  by  four.   The  plotter  and 
paper  should  be  readied  before  the  keyboard  input  data  are  read. 


105 


FORTX 

Synopsis 

This  1800  program  computes  the  Fourier  Transform  of  the  thermal 
neutron  pulse  stored  on  disk  by  TKERM.   The  integration  scheme  used  is 
Simpson's  Rule,  with  the  initial  channel  contribution  being  computed 
by  the  trapezoidal  rule.   The  sumraation  increment  is  greater  than  the 
channel  width  by  an  amount  equal  to  the  deadtime  betv;een  channels. 
The  midpoint  of  the  summation  increment  coincides  with  the  midpoint 
of  the  channel.   The  frequency  range  and  interval  arc  specified  in  the 
input  data.   The  frequency,  amplitude  and  phase  are  punched  on  paper 
tape  and  printed. 

Input 

The  following  data  are  read  from  a  card  in  FORMAT  (3I-,  SSlO.ij-): 
LMIN  =  file  number  of  first  pulse 
LMAX  =  file  number  of  last  pulse 
NCHAN  =  nunber  of  channels 

TMCCH  =  multichannel  analyzer  channel  width  in  microseconds 
TMCDT  =  deadtime  between  channels  in  microseconds 
FRiMIN  =  initial  frequency  (cps)  of  Fourier  analysis 
FRMAX  =  final  frequency  (cps)  of  Fourier  analysis 
DELTA.  =  frequency  interval  (cps)  . 

Output 

The   printed   output    includes,    at  each   frequency,    the   real   and 
im.aginary  components   of   the    transformed  pulse,    the   amplitude   and   the 
phase.      The   paper    tape   output   includes   only   the   frequency,    amplitude 
and  phase. 
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Special  Controls 

PAUSE  ^^3:5.--   This  pause  occurs  when  the  analysis  of  the  specified 
pulses  has  been  completed.   Pressing  START  returns  the  program  to  the 
reading  of  the  data  input  card. 

Comnents 

The  paper  tape  output  is  punched  for  input  into  the  ALPllA  and  XI 
Codes  which  will  be  described  next.  One  frequency  analysis  requires 
approximately  20  seconds  when  1,021+  channels  of  data  are  used. 

ALPHA 

Synopsis 

This  l80O  prograr.i  plots  the  experimental  values  of  the  amplitudes 
as  a  function  of  distance  for  each  frequency  on  a  semilog  scale.   A 
rough  plot  for  even  frequencies  only  is  given  so  that  a  check  can  be 
made  for  asy^r.ptotic    decay  and  bad  data  points.   A  second  plot  is  then 
made  which  plots  the  dat?.  points  for  all  frequencies,  computes  a 
semilog  least  squares  fit  of  the  form 

F(z)  =  ae^^ 

where  a  and  b  are  fitting  parameters  and  F(2)  is  the  experimental 
data  and  draws  the  least  squares  fit  through  the  data  points.   The  ex- 
perimental points  to  be  included  in  the  connputation  are  fixed  from 
1  to  M-VX  or  can  be  varied  by  using  the  sense  switches  and  data  input. 

Input 

The  output  tape  cf  F01\.TX  is  required  input.   51  frequency  data 
are  read  for  each  data  position  (0-500  cpSj  in  steps  of  10  cos).   The 
following  data  are  read  from  the  keyboard  in  FOS>!AT  (Fli.O,  l6li;)  : 
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YMIN  =  minimum  value  to  be  plotted  (should  be  an  integer 

povar  of  10) 
NCYCL  =  number  of  cycles  to  be  used  along  vertical  axis 
IMIN  =  1 
IMX  =  number  of  data  runs  to  be  included  in  the  least 

squares  fitting 
MAXFR  =  one  plus  the  number  of  frequencies  to  be  analysed 

(frequencies  are  incremented  by  10  cps) 
ILLFR  =  the  number  of  the  frequency  at  which  data  elimination 

is  to  begin 
ISTEP  =  minimum  number  of  data  points  to  be  considered  in  the 

progressive  scanning  of  the  data 
IBAD  =  data  run  numbers  to  be  omitted,  maximum  of  10,  leave 

blank  if  none. 

Output 

If  the  data  are  being  scanned  the  printed  output  includes  the 
initial  and  final  data  positions,  the  asymptotic  decay  constant,  the 
zero  value  and  the  deviation  of  each  from  the  previous  computation. 
Any  om.itted  data  poistions  are  also  listed.   If  the  data  are  not  being 
scanned,  the  deviations  are  not  printed. 

UTien  there  is  no  scanning  and  the  ILLF^lth  frequency  is  about  to 
be  computed  a  message  is  printed  to  input  the  initial,  final  and 
intermediate  data  positions  to  be  omitted  at  the  keyboard.   This  pro- 
vision allows  the  omission  of  any  run  from  the  computations. 

The  m.inimum  number  of  data  runs  allo'.;ed  is  three.  If  less  than 
three  are  being  considered,  che  computation  is  omitted  and  the  next 
frequency  cor.putation  is  .T.ade.   A  "NO  A3Y?-ffTCTIC  DECAY"  and  frequency 


103 

identification  ir.essage  is  printed  in  this  case. 

The  frequency  and  syniptotic  decay  constants  are  punched  on  paper 
tape  at  the  end  of  the  program.   This  tape  is  an  input  to  the  PL0T2 
program  which  will  be  described. 
Special  Controls 

SENSE  SWITCH  0.--If  SENSE  SWITCH  0  is  turned  on,  no  plots  are 
ma  de . 

SENSE  SWITCH  l.--If  SENSE  SWITCH  1  is  turned  on,  no  printed  out- 
put is  obtained. 

SENSE  SWITCH  3. --If  SENSE  SWITCH  3  is  on  the  data  points  are 
scanned.   The  scanning  procedure  is  such  that  ct  is  computed  and 
printed  for  the  first  ISTEP+1  data  points.   The  group  of  ISTEP+1 
points  is  then  shifted  over  one  data  position  and  a  is  computed  again. 
These  processes,  computing,  printing  and  shifting,  are  continued 
until  the  last  ISTSP+1  data  points  are  considered.   This  completes  one 
scan  of  the  data.   The  next  scan  is  made  with  ISTEP  increased  bv  one. 
Additional  scans  are  made  until  the  group  of  data  points  considered 
equals  the  total  number  of  data  points 

PAUSE  2i2;f.--This  pause  occurs  so  that  the  paper  on  the  plotter 
may  be  changed  between  the  rough  and  final  plots. 

PAUSE  1111. --This  pause  occurs  at  the  completion  of  the  cor.iputa- 
tions.  Pressing  START  returns  the  prograr.i  to  the  initial  reading  of 
the  keyboard  input  data. 

Comments 

Several  ruas  of  ALPH.^  may  be  necessary  to  obtain  the  proper 
asymptotic  region  and  to  determine  which  data  points  are  to  be  elim- 
inated.  For  these  runs  the  tape  output  can  be  suppressed  by  turning 
the  paper  tape  punch  off. 
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XI 
Synopsis 

This  1800  program  is  very  similar  in  logic  to  the  ALPHA  program. 
However,  a  linear  least  squares  fit  is  made  to  the  corrected  phase 
shifts  as  a  function  of  frequency.   The  correction  is  required  be- 
cause the  range  of  arguments  for  inverse  angles  is  from  -  /2  to  +  /2. 
Arguments  less  than  -  /2  must  be  corrected  by  subtracting  the 
appropriate  multiple  of  n .   This  correction  is  performed  by  XI.   A 
tape  of  the  frequencies  and  corrected  phases  can  be  obtained. 

Input 

The  data  output  tape  of  FORTX  is  also  read  by  XI.   The  following 
numbers  are  read  from  the  keyboard  in  FOR>Lf\T  (Fl+.O,  l^X^)  : 

YWAX  =  the  maximum  phase  shift  (lag)  to  be  plotted 

IMIN  =  1 

IMAX  =  the  number  of  data  runs 

MAXFR  =  one  plus  the  number  of  frequencies  to  be  analyzed 

(frequer.cies  are  incremented  by  10  cos) 
ILLFR.  =  the  first  frequency  number  at  vhich  data  elimination 

is  to  begin 
ISTEP  =  initial  number  of  points  considered  in  the  scanning 

of  the  data  positions 
IBAD  =  the  data  runs  to  be  omitted,  maximum  of  10,  leave 

blank  if  none. 

Output 

The  outnut  of  XI  is  identical  to  the  output  of  ALPfL\. 
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Special  Controls 

SENSE  SWITCH  2 .--If  SENSE  SWITCH  2  is  turned  on,  the  corrected 
phases  are  punched  on  paper  tape. 

Comments 

The  comments  for  ALPHA  are  applicable.   Generally  ALPHA  and  XI 
are  run  successively.   If  several  runs  of  XI  are  required,  the  paper 
tape  punch  can  be  turned  off  until  the  final  run.   Then  a  continuous 
tape  input  to  PL0T2,  to  be  described  next,  is  obtained. 

PL0T2 

Synopsis 

PL0T2  is  an  l800  program  which  plots  a  versus  frequency  (cps)  ,  | 

2 
versus  frequency  (cps),  |  versus  a,  2Cis  versus  a>  (rad/sec)  ,  and  oc     - 

2         2  2   2 

I  versus  co  .   20:5  and  ^  -  |   are  computed  and  printed  as  well  as 

plotted. 

Input 

The  ALP-XI  tape  which  is  the  ALPHA  and  XI  Codes  output  is  read  at 
the  tape  reader.   The  following  data  are  read  fro~  card  in  FORi-'iAT 
(2F4.0,  51^): 

X>LAX  =  maximum  value  of  horizontal  axis 

YMAX  =  maximum  value  of  vertical  axis 

MAXFR  =  maximum  number  of  frequencies 

NO-    =  point  plot  s^,njibol  control  and  run  number  for  UFNLLS 


ICARD  =  card  punch  option;  if  0,  no  cards  punched;  if  1,  2 


2    2 
and  OL     -    c,      are  punched  out  for  input  to  UFNLLS  Code, 


Ill 


Output 

2    2 
The  value  of  2ai   and  a  -  ^   are  printed.   Depending  on  the  value 

of  ICAR.D,  cards  are  punched  for  UFNLLS  input.   The  independent  variable 

2    2  2 

on  the  a  -  ^  cards  is  od  rather  than  od.   The  five  plots  mentioned 

in  the  synopsis  are  also  obtained. 

Special  Controls 

PAUSE  1111. --This  pause  occurs  so  that  the  ALP  tape  output  of 
ALPHA  can  be  readied  at  the  tape  reader.   When  readied,  press  START 
to  read  the  tape. 

PAUSE  2222. --This  pause  occurs  after  the  reading  of  the  ALP  tape. 
It  occurs  so  that  the  XI  tape  output  can  be  readied  at  the  tape  reader. 
When  readied,  press  START  to  read  the  tape. 

PAUSE  2222. --These  pauses  occur  between  the  plots  to  allow  the 
paper  on  the  plotter  to  be  readied  for  the  succeeding  plots.  When 
readied,  press  STAXT  to  continue  the  plots. 

PAUSE  33^3. --This  pause  occurs  when  all  plots  have  been  made. 
Pressing  START  returns  the  progran  to  the  reading  of  the  data  input 
card. 

ConTrnents 

The  maximuni  values  read  in  for  the  vertical  and  horizontal  axes 

apply  only  to  the  U   and  5  axes.   The  maxi:nuri  on  the  frequency  axis  is 

2    2 
14.00  cps.   The  maxima  on  the  vertical  axes  for  2a|  and  0!  -  =   are 

2  5 

0.02  and  0.0^5  respectively,   cu  and  00  maxima  are  set  at  2.^x10  and 

5.0x10  respectively.   Should  it  be  dasired  to  change  any  of  th-ese 

latter  values,  the  appropriate  scaling  factors  raust  be  changed  and  the 

progran  recor.piled. 
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PLOT^ 

Synopsis 

This  1800  program  plots  and  prints  the  theoretical  values  of  a 
and  ^.   The  FERVE2  Code  punches  the  value  of  frequency,  a  and  |  on 
cards  for  input  to  PL0T3.   The  theoretical  values  can  be  point  or 
line  plotted.   Options  are  available  for  obtaining  theoretical  plots 
which  are  drawn  over  the  experimental  point  plots  obtained  from  the 
PL0T2  Code. 

Input 

The  following  data  arc  read  from  card  in  FORi-L-.T  (,^lk)  : 
WiX     =   nur^ber  of  frequencies 
NO   =  point  plot  symbol  control 

NRUN  =  axis  drawing  and  labeling  control;  if  0  axes  are  drawn 
and  labeled,  if  not  0  only  the  theoretical  values  are 
plotted 
LORPT  =  line  or  point  plot  control;  1  =  line,  2  =  point  plots 
NDUN  =  plot  control;  if  =  1,     |  versus  ci   is  plotted;  if  -  2, 
CC  versus  frequency,  '-   versus  frequency  as  well  as  5 
versus  a:  are  plotted. 
Frequency,  y.   and  :  are  read  frcn  card  in  FCRI^'AT  (5EI2.5).   MAX  cards 
are  read. 

Output  . 

In  addition  to  the  graphs  the  theoretical  values  of  a   and  |  are 
printed,  depending  on  SENSE  SWITCH  1. 
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Special  Control 

PAUSE  1111. --This  pause  occurs  so  that  the  plotter  may  be  readied 
for  the  a  versus  frequency  plot.   Push  START  when  readied  to  obtain 
plot. 

PAUSE  2222. --Same  as  PAUSE  1111  except  that  ^  versus  frequency 
is  plotted. 

PAUSE  ^5^^. --Same  as  Pause  1111  except  that  |  versus  a  is  plotted. 

PAUSE  l4.1j.lj.l4..--xhis  pause  occurs  so  that  the  printer  may  be  readied. 
When  readied  press  ST.'\RT  to  continue. 

PAUSE_^2^. --Press  START  to  continue. 

SENSE  SWITCH  l.--I£  SENSE  SV/ITCH  1  is  turned  on  the  listing  of 
theoretical  values  of  a  and  |  will  not  be  made. 

Coniraents 

If  only  the  listing  is  desired  it  can  be  obtained  by  setting 
NDUN  =  1  and  turning  SENSE  SWITCH  0  on. 

The  axes  of  the  a  versus  frequency  and  |  versus  frequency  plots 
can  not  bs  drawn  or  labeled  by  PLOT^.   These  are  drawn  by  PL0T2  and 
are  positioned  on  the  plotter  so  that  the  theoretical  plots  are  made 
over  the  experimental  plots.   The  plots  are  positioaed  by  taping  the 
plots  on  the  drum  and  placing  the  alignment  reticle  over  the  origin. 
The  san-.e  holds  true  for  the  j  versus  a  plots  in  that  the  reticle  n;ust 
be  placed  over  the  origin  of  the  PL0T2  plot  if  the  plots  are  to  be 
superimposed.   The  PL0T2  plots  and  PLOTJ  plots  are  o_ffse_t  by  one  inch 
in  the  horizontal  direction. 

FERVE2 

Synopsis 

This  560/50  program  computes  the  theoretical  values  of  a.   and  \ 


from  equation   [2.12]   with  n,p   =   1   and  2.      The   Newton-Raphson   itera- 
tion  technique    is   us.d   to   obtain   solutions.      At   each   frequency,    initial 
estimates   for  a  and   |    are   read. 

Input 

The  first  card  read  is  a  title  card,  FORMAT  (12A6) .   Three 
initial  data  input  cards  are  read  for  each  set  of  computations.   The 
following  data  are  read  in  F0R2-LVT  (6E12.6): 

BTSQl  =  fundamental  transverse  buckling  squared 

BTSQ2  =  first  higher  transverse  buckling  squared 

SQLO  =  diffusion  length  squared 

TAU   =  age  to  thermal 

DTH   =  diffusion  constant 

VTH   =  thermal  velocity 

SL    =  slowing  down  time 

AO    =  alpha  zero 

CI    =  cooling  index 

GAMMA  =  Feinberg-Galanin  constant 

PESCF  =  resonance  escape  probability 

ETA   =  fast  neutrons  produced  per  thermal  neutron  absorption 

JONE   =  J  (B  *RTOT)  where  B,  =  BTSQl 
o^  1     -^        1      ^ 

JTVJO  =-  J  (B  *RTOT)  where  B,  =  BTSQ2 
RTOT  =  R,  extrapolated  radius  of  assembly 
RZERO  =  radius  of  fuel  rod 
Next,  the  following  data  are  read  in  FORMAT  (I1|,E12.6): 
MAXIT  =  maximum  number  of  interations 
-  E     =  tolerance  of  computing  a  and  P 
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The  last  set  of  cards  read  contain  the  following  date  in  FORl-lAT 
(5F12.6,  211+)  : 

03     =  frequency  (cps) 

ALPHA  =  initial  estimate  of  a 

PSI   =  initial  estimate  of  | 

ILAST  =  0,  for  last  card  in  the  set 

=  1,  for  all  other  cards  in  the  set 

ISENT  =  0,  for  last  card  in  the  last  set 
=  1,  for  all  other  cards 
Succeeding  sets  must  have  each  of  the  above  cards  included. 

Output 

The  title  card,  initial  input  data  and  the  frequency  analysis  at 
each  frequency  is  printed  out.   Cards  of  frequency,  a  and  '  are 
punched  for  input  to  PLOTJ. 


APPENDIX    E 
EXPEPvIMENTAL    AMPLITUDES    AND    PH^^.SES 
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HEAVY   WATER 
Z    =    50.     CflS    AXIAL    POSITION 


FRrQ(CPS) 

AMPLITUDE 

PHASE (RADIANS) 

0. 

0.207609F. 

03 

O.OOOOOOE  00 

10. 

0.20it8U3E 

05 

0.335188E  00 

20. 

0.197036E 

03 

0.653727E  00 

30. 

0.185625E 

03 

0.9  8  0767E  00 

uo. 

0.171963E 

05 

0.12837'tE  01 

50. 

0.157391E 

03 

-0.1569ii9E  01 

60. 

0.142318E 

03 

-0.129513E  01 

70. 

0.128813E 

03 

-0.103356E  01 

80. 

0.115573E 

03 

-0.7S3382E  00 

90. 

0.1035U1E 

03 

-0.5U33'j2E  00 

100. 

0.92Uti02E 

02 

-0.312389E  00 

110. 

0.8235S1+E 

02 

-0.896399E-01 

120. 

0.732505E 

02 

0.12570!+E  00 

130. 

0.650U37E 

02 

0.33U251E  00 

140. 

0.57691SE 

02 

0.535i;0t+E  00 

150. 

0.511135E 

02 

0.732852E  00 

160. 

O.Jt523U5E 

02 

0.9233SiiE  00 

170. 

0.lt0n039E 

02 

0.1109C5E  01 

130. 

0.353537E 

02 

0.1291'4'tE  01 

100. 

0.312251E 

02 

0.1U5912E  CI 

200. 

0  .  2  7  5 1+  J  9  E 

02 

-0.1I>97S3E  01 

210. 

0.2U2528E 

02 

-0.1326U0E  01 

220. 

0.215231E 

02 

-0.115S32E  01 

230. 

0.186967E 

02 

-0.99tiOQ6E  00 

21+0. 

0.163613E 

02 

-J.S3tt229E  00 

250. 

0.1I+3035E 

02 

-0.5790ti3E  00 

260. 

0  .  1 2  5 1 5  3  E 

02 

-0.527121JE  00 

270. 

0.109[+5!+E 

02 

-0.376865E  00 

2S0. 

0.955673E 

01 

-0.227303E  00 

290. 

0.832303E 

01 

-0.799290E-01 

300. 

0.722U35E 

01 

0.555321E-01 

310. 

0.625999E 

01 

0.203637E  00 

320. 

0.5U1C12E 

01 

0.3512iilE  0  0 

330. 

0.lt657ii2E 

01 

0.  ft 95317 E  00 

3!+0. 

0.395737E 

01 

0.55502!+E  0  0 

550. 

0  .  3  3 '.)  1  3  3  E 

01 

0.775390E  00 

lis 


HEAVY  WATER 
=    55.    CMS    AXIAL    POSITION 


REQ(CPS) 

AMPLITUDE 

PIIASE(RADlAf!S) 

0. 

0.172123E 

03 

0.000900E 

00 

10. 

0.159659E 

03 

U.355192E 

00 

20. 

0.162780E 

03 

n.705325E 

00 

30. 

0.152622E 

03 

0.103315E 

01 

itO. 

0.1ii0552E 

03 

0.155987E 

01 

50. 

0.12777t+E 

05 

-0.1I»7631E 

01 

GO. 

0.11506SE 

0  5 

-0.113701E 

01 

70. 

0.1029U1E 

05 

-0.911175E 

0  0 

80. 

0.916b95E 

02 

-0.6t*7927E 

00 

90. 

0.813702E 

02 

-0.395321E 

00 

100. 

0.720602E 

02 

-0.153i+fi3E 

00 

110. 

.     0.63G926E 

02 

0  .  8  0  5  9  2  9  E  - 

-01 

120. 

0.5S1970E 

02 

0.305521E 

00 

150. 

0.I49U9U8E 

02 

0.5259ttUE 

00 

lUO. 

0.I+55153E 

02 

0.73S7U1E 

00 

150. 

0.581S31E 

02 

0.9i451!jUE 

0  0 

160. 

0  .  5  3  It  7  3  S  E 

02 

0.11U550E 

01 

170. 

0  .  2  9  3 1 6  0  E 

02 

0.135997E 

CI 

180. 

0  .  2  5  G  5 1+  5  E 

02 

0.152960E 

01 

190. 

0  .  2  2  ij  3  2  1 E 

02 

-  0  . 1 1+  2  G  5  5  E 

01 

200. 

0.195S53E 

02 

-0.12!i511E 

01 

210. 

0  . 1  7  0  7  3 1 E 

02 

-0.106765E 

01 

220. 

0.1f+3553E 

02 

-O.S3'+G70E 

00 

250. 

0.1291 72 E 

02 

-0.726173E 

0  0 

21+0, 

0  . 1 1 2  2  U  3  £ 

02 

-0.561105E 

GO 

250. 

0.974399E 

01 

-0. 39393 3E 

0  0 

260. 

0.8U5716E 

01 

-0.233't21E 

00 

270. 

0.75203itE 

01 

-0.802301+E- 

-01 

2  80. 

0  .  6  3  2  3  6  7  E 

01 

0.763991E- 

-01 

2  00. 

Q.5'-'rhh2klL 

01 

0.230392E 

CO 

500. 

0.1467342E' 

01 

0.58325bE 

00 

510. 

0.399 330 E 

01 

0.535261E 

00 

520. 

0.3ia302E 

01 

0  .  6  8  1  5  !i  7  E 

00 

550. 

0  .  2  3  9  9  9  7  £ 

01 

0  .  8  3  0  0  U  2  E 

00 

3U0. 

0.2'jU303E 

01 

0.976290E 

00 

350. 

U  .  2  0  5  5  5  3  F 

01 

0  .  1 1 1 9  5  S  E 

01 
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HEAVY  WATER 
Z    =    50.    CMS    AXIAL    POSITION 


FREQ(CPS) 

AilPLITUDE 

PHASE (RADIANS) 

0. 

0.140965E 

03 

O.OOOOOOE 

00 

10. 

0.133 8 30E 

0  3 

0.3757ii8E 

00 

20. 

G.132839E 

03 

0.7itH057E 

00 

30. 

0  .  1 2  fj  0 1  9  E 

03 

0.109919E 

01 

hO. 

0.113592E 

03 

0  . 1  !t  3  3  0  9  E 

01 

50. 

0  .  1  0  2  5  9  0  E 

03 

..     -0.133167E 

01 

60. 

0.917itii9E 

02 

-0.1076I(1E 

01 

70. 

C.81U99iiE 

02 

-0.786362E 

00 

80. 

0.720714E 

02 

-0.509335E 

00 

90. 

0.635330E 

02 

-0.2U522iiE 

00 

100. 

0  .  5  5  3  7  5  2  E 

02 

0.887002E- 

-02 

]10. 

0.ft90492E 

02 

0.25366SE 

00 

120. 

0  .  U  2  9  8  3  3  E 

02 

0.i<8995!tE 

00 

130. 

0.376205E 

02 

0.718389E 

00 

140. 

0  .  3  2  3  8  6  7  E 

02 

n.939G09E 

00 

150. 

0.2872Q0E 

02 

0.115't02E 

01 

160. 

0  .  2  5  0  3  0  3  E 

02 

0  .  1  3  6  2  7  5  E 

01 

170. 

0.213 74 5E 

02 

0.155509E 

01 

130. 

0.ig057GE 

02 

-0.137G99E 

01 

190. 

0.165751E 

02 

-0.113  3UUE 

01 

200. 

0.U3959E 

02 

-0.9950[|2E 

00 

210. 

0.12US55E 

02 

-0.S1135UE 

00 

220. 

0.10323i+E 

02 

-0.631H+05E 

00 

230. 

0.933I+93E 

01 

-0.I+61325E 

00 

2tj0. 

0.813592E 

01 

-0.29U37E 

00 

250. 

0.70ttl55E 

01 

-0.125053E 

00 

260. 

0.50373UE 

01 

0.363151E- 

■01 

270. 

0.52G0S7F 

01 

0.194275E 

00 

230. 

0.I455O35E 

01 

0.3'+7203E 

0  0 

290. 

0.39U502E 

.01 

0.f+9760GE 

00 

300. 

0.3U2835E 

01 

0.6I+3100E 

00 

310. 

0.293131E 

01 

0.801975E 

00 

320. 

0.2533U7E 

01 

0.960259E 

00 

330. 

0.223521E 

01 

0.112698E 

01 

3i+0. 

0  .  1 9  0  2  7  0  E 

01 

0.13022SE 

01 

350. 

0  . 1 5  3  5  5  U  E 

01 

0.U7523E 

01 

120 


Z    = 


FRF.Q(CPS 

0. 

10. 

20. 

30. 

liO. 

50. 

GO. 

70. 

80. 

90. 

100. 

110. 

120. 

130. 

lUO. 

150. 

160. 

170, 

130. 

190. 

2  00. 

210. 

220. 

230. 

2U0. 

250. 

25  0. 

270. 

2  SO. 

290. 

300. 

310. 

320. 

330. 

3U0. 

350. 

HEAVY  WATER 

55.     CMS    /! 

'.X!AL 

POSITION 

AMPLITUDE 

PHASE(RADIAMS) 

0.1102257. 

03 

O.OOOOOOE 

00 

0.11U360E 

03 

0.395917E 

00 

0.10915itE 

03 

0.784036E 

00 

O.lOlUgiiE 

03 

0.115820E 

01 

0  .  9  2  U  7  7  0  E 

02 

0.151509E 

01 

0.83015SE 

02 

.     -0.128777E 

01 

0.737502E 

02 

-0.955G914E 

00 

0.6505b9E 

02 

-0.6G175GE 

00 

0.5710't5E 

02 

-0.371212E 

00 

o.ttggsoor: 

02 

-0.937075E- 

■01 

0.U359f*7E 

02 

0.17225SE 

00 

0.379825E 

02 

0.ft23135E 

CO 

0.330U55E 

02 

0.67I4SU2E 

00 

0.287173E 

02 

0.913193E 

00 

0  .  2  (+  9  3  2  2  E 

02 

0.11it37SE 

01 

0.2163^9E 

02 

0.13G715E 

01 

0.137703E 

02 

-0.155729E 

01 

0.162792E 

02 

-0.13U553E 

01 

0.1I+1105E 

02 

-0.113352E 

01 

0.1221I+5E 

02 

-0.935375E 

00 

0  . 1 0  5  6  0  5  E 

02 

-0.737652E 

00 

0.911U63E 

01 

-0.5I+3U25E 

00 

0.785329E 

01 

-0.353527E 

00 

0.67503:iE 

Gl 

-0.1672U0E 

00 

0.573233E 

01 

0.13551+5E- 

-31 

O.U9[+559E 

01 

0.1S71!;6E 

00 

0.ii23b32E 

01 

0.35356SE 

00 

0.36U[;'aE 

01 

0  .  5  1 5  9  7  9  E 

00 

0.314741E 

01 

0.673233E 

0  0 

0.2715t;3E 

01 

0.811G652F 

00 

0.233118E 

01 

0.101526E 

01 

0.19  9  019E 

01 

0.113ii35E 

01 

0  .  1 6  3  6 14  5  E 

01 

0.13I43C7E 

01 

0.1U3023E 

01 

0.150313E 

01 

0  . 1 2  1 0 1  3  E 

01 

-0.1U7371E 

01 

0.1025G0E 

01 

-0.131550E 

01 
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HEAVY   WATER 
Z    =    70.    C!1S    AXIAL    POSITIOri 


FREQ(CPS) 

AMPLITUDE 

PHASE (RAO lANS) 

0. 

0.97ia31F. 

02 

O.OOOOOOE 

00 

10. 

0.95773ir: 

02 

0.'4l6G27E 

00 

20. 

0.911SliiE 

02 

0 . 8 2 U 9  CI E 

00 

30. 

Q.ZhkSSkE 

02 

0.12187iiE 

01 

1^0. 

0.7G5G;i9E 

02 

-0.15i+729E 

01 

50. 

0.683255E 

02 

-0.119108E 

01 

60. 

0.603079E 

02 

-0.S53507E 

00 

70. 

0.523364E 

02 

-0.533320E 

00 

80. 

0.460533E 

02 

-0.22  8!i55E 

00 

90. 

O.ftOOlOlE 

02 

0.527SS8F- 

-01 

100. 

0.3ttG650E 

02 

0  .  3  U 1  3  9  0  E 

00 

110. 

.  0.299733E 

02 

0.609923E 

00 

120. 

0.253783E 

02 

0.8G774GE 

00 

130. 

0.223220E 

02 

0  . 1 1 1 6  2  2  E 

01 

11*0. 

0.192U00E 

02 

0.135G16E 

01 

150. 

0.16577GE 

02 

-0.155372E 

01 

IGO. 

0  . 1 U  2  9  2  5  E 

02 

-0.132950E 

01 

170. 

0.12334^E 

02 

-0.11I99SE 

01 

130. 

0.10Gt*2iiE 

02 

-  0  .  S  9  6  9  3  7  E 

00 

190. 

0.917352E 

01 

-0.687352E 

00 

200. 

0.7892:+0E 

01 

-0.U32U26E 

00 

210. 

0.6 779  GIF 

01 

-0.2C3570E 

00 

220. 

0.5S2571E 

01 

-0.9054U5E- 

-01 

230. 

0.50n935E 

01 

0.930G35E- 

-01 

2ii0. 

0  .  U  3  0  3  7  9  E 

01 

0.233727E 

00 

25U. 

0.3707U3E 

01 

0.it673U5E 

00 

2u0. 

0.313U97E 

01 

0  .  6  5  0  2  3  U  E 

0  0 

270. 

0.272UU2E 

01 

0.S31759E 

00 

2  ;>  0 . 

0  .  2  3 1 S  S  U  E 

01 

0.1C09CGE 

01 

29  0. 

0.19GC02E 

01 

0.1 179  0  IE 

01 

300. 

0.167U2:E 

01 

0.13ii309E 

01 

310. 

0.1I42S61E 

01 

0.150!+72E 

01 

320. 

O  .  1 2  2  C  0  1  E 

01 

-0.1t47232E 

01 

^  7  <-% 

J  5u  , 

0.10571ttE 

01 

-0.12931UE 

01 

3U0. 

0  .  3  9  3 1 2  1 E 

00 

-0.109U91E 

01 

350. 

0  .  7  3  9  7  S  3  E 

00 

-  0  .  S  2  3  9  2  3  E 

00 

122 


HEAVY  WATER 
Z    =    75.    CMS    AXIAL    P05ITI0M 


FREQ(CPG) 

AMPLITUDE 

PHASE(RAOIAr;S 

0. 

0.812i!lOE    02 

O.OOOOOOE    00 

10. 

0.7980i»3E    02 

0.435019E    00 

20. 

0.757911F    02 

0.861!ii*2E    00 

30. 

0.C99291E    02 

0.127263E    01 

hO. 

0.5307itOE    02 

-0.1Ii76dSE    01 

50. 

0.559«i02E    02 

-0.110i+l48E    01 

60. 

0.U90253E    U2 

-0.751889E    00 

70. 

0.'*2G1G3E    0  2 

-0.m7502E    00 

80. 

0.368f401F    02 

-0.995358E-01 

90. 

0.31729t»E    02 

0.20379GE    00 

100. 

0. 2 72 5 72 E    0  2 

O.U914222E    00 

no. 

0.233701E    02 

0.77311ftE    00 

120. 

0.200071E    02 

O.lQfflSSE    01 

130. 

0.171079E    02 

0.130032E    01 

lUO. 

0.1f»G135E    02 

0.155032E    01 

150. 

0. 12 1; 819 E    02 

-0.13ii92  6E    01 

IGO. 

0.10bU97E    02 

-0.111i;3  2E    01 

170. 

0.90o02  5E    01 

-0.333017E    GO 

ISO. 

0.77U73oE    01 

-0. 663 3 95 E    00 

190. 

0.5G2309E    01 

-0.ft55'^30E    00 

200. 

0.56G95oE    01 

-0.2'f5ft70E    0  0 

210. 

0.US5076E    01 

-0.372373E-01 

220. 

0.iil3533E    01 

0.163 79 3E    0  0 

230. 

0.350710E    01 

0.369ft';lE    00 

2f+0. 

0.29671i+F    01 

0. 551 59 3 E    GO 

250. 

0.251293E    01 

0.7itG155E    00 

2G0. 

0.213242E    01 

0.9253i+3E    00 

270. 

0.131I+5I+E    01 

n.ll0296E    01 

280. 

0.15  395SE    01 

n.l2799SE    01 

290. 

0  .  1 3  li  3  1  7  E    01 

0.i£t52it0E    01 

300. 

O.IIOIOIE    01 

-0.152213E    01 

310. 

0.93it959E    00 

-0.136139E    01 

320. 

0.797i52E    00 

-0.119!+55E    01 

330. 

0.675799E    00 

-0.102I+36E    01 

5k0. 

0.570312F    00 

-0.352211E    00 

3  5  0 . 

0.U7323UE    00 

-0.6  735G3E    0  0 
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HEAVY   WATER 
30.    CMS    AXIAL    POGITION 


FRF.Q(CPS) 

AMPLITUDE 

PIIASFCRADI/' 

^MS) 

0. 

0.671335E 

02 

O.COOOCOE 

00 

10. 

0.659071F 

02 

0.ti53303F 

00 

20. 

0.62U739E 

02 

0.897772E 

00 

30. 

0.57l|G37r: 

02 

U.132655F 

01 

liO. 

0.51G2tjOE 

02 

-0.1'»0579E 

01 

50. 

0.ti55Gl'*n 

02 

-0.1017f+i;E 

01 

60. 

0.397096E 

02 

-0.5U9577F 

00 

70. 

0.3U310I)E 

02 

-0.50n7G5E 

00 

80. 

0.2Q'^729F 

0  2 

0.5073ottE- 

-01 

90. 

0  .  2  5  2  1 5  8  F 

02 

0.34G922F 

00 

100. 

0  .  2  1 5 1 2  2  F 

02 

0.6it9199E 

00 

110. 

.  0.183212E 

02 

0.933370E 

00 

120. 

0.155378F 

02 

0.121755F 

01 

130. 

0.132505F 

02 

0.1U8550F 

01 

lUO. 

0  . 1 1 2  5  1  3  F 

02 

-0.139fi85E 

01 

150. 

0.95^000E 

01 

-0.11!4352F 

01 

IGO. 

0.809025F 

91 

-0.900912F 

00 

170. 

0.58G228F 

01 

-0.6o5555F 

00 

180. 

0.582559E 

01 

-  0  .  4  3  G  7  3  G  E 

00 

I'JO. 

0.l*9i+20  0E 

01 

-0.212^3GF 

00 

200. 

O.itlOOOOF 

01 

O.SOOttSlE- 

-C2 

210. 

0.35501SF 

01 

0.220217E 

00 

220. 

0.30021GE 

01 

0.f+29533F 

00 

230. 

0.253Utt4E 

01 

0.5351GUF 

00 

2U0. 

0.213005E 

31 

0  .  S  3  5  4  5  It  E 

09 

250. 

0  .  1 7  C  i»  U  3  E 

01 

0  . 1  0  2  3  3  G  E 

01 

2G0. 

0.1i+3757E 

91 

0.121273F 

01 

270. 

0.12i;276F 

01 

0.1532lt3F 

01 

220. 

0.10i|GU2E 

01 

0  .  1 5  4  2  9  7  F 

91 

290. 

0.890191E 

0  0 

-  0  .  1 1;  1+  9  1  7  F 

01 

500. 

0.759523E 

0  0 

-0.127905F 

01 

310. 

0.650979F 

00 

-0.111S06F 

01 

320. 

0.558G2itE 

00 

-0.9tt3535E 

00 

350. 

0.I+75525F 

90 

-0.7733i+2F 

on 

3U0. 

0.399tj5SF 

0  0 

-  0  .  5  0 1 3  G  2  F 

0  0 

350. 

0  .  3  3  5 1  7  S  E 

00 

-O.i+f+593'tF 

00 

12l| 


HEAVY   WATER 
Z    =    85.    CilS    AXIAL    POSITION 


fr[:q(cps) 

AMPLITUDE 

PHASE(RADIAf!S) 

0. 

0.5'tit303E 

02 

O.OOOOOOE 

CO 

10. 

0.533918E 

02 

0.ii7177'tE 

00 

20. 

0.50U930E 

02 

0.93ti35CE 

00 

30. 

0.i*G276iiE 

02 

0.13S067E 

01 

kO. 

0  .  tt  1 3  G  S  9  E 

02 

-0.133'»96E 

01 

50. 

0.3G3021E 

02 

-0.930807E 

00 

GO. 

0  .  3 1  ii  3  7  5  E 

02 

-0.5'i3011E 

00 

70. 

0.2G9763E 

02 

-0.185107E 

00 

30. 

0.25001GE 

02 

0.159G99E 

CO 

90. 

G.19523tiE 

02 

O.ftSSlGSE 

00 

100. 

0.165320E 

02 

O.S020GOE 

00 

110. 

0.139G7GE 

02 

0.110271E 

01 

120. 

0.117SS5E 

02 

0  . 1 3  9 1 G  2  E 

01 

130. 

0.99'iOOU: 

01 

-0.1if7167E 

01 

lUO. 

0.837901E 

01 

-0.120301E 

01 

150. 

0.705377E 

01 

-0.9tt23U2E 

CO 

IGO. 

0.59U516E 

01 

-0.5911f+3E 

00 

170. 

0.5007it9E 

01 

-0.Utt6332E 

00 

180. 

0  .  4  2  1 3  5  9  E 

01 

-0.2087!+'*E 

00 

190. 

0  .  3  5  U  G  2 1;  E 

01 

0.21ftlGSE- 

-01 

200. 

G.29  33  0  3E 

01 

0  .  2 ';  ;j  G  2  9  E 

00 

210. 

0.25253iiE 

01 

0  .  I*  G  2  3  9  1 E 

00 

220. 

0.213720E 

01 

0.679G57E 

00 

230. 

0.180753E 

01 

0.894609E 

0  0 

21*0. 

0.1526U0E 

01 

O.llOOOGE 

01 

250. 

0.123280E 

01 

0  .  1 3  2  1 5  2  E 

01 

260. 

0  .  1 0  7  5  3  8  E 

Ul 

0.153009E 

01 

270. 

0.897U15E 

0  0 

-0.140352: 

01 

280. 

0. 743  0 30 E 

OD 

-O.12OI714E 

01 

290. 

0.612852E 

0  0 

-0.10092:iE 

01 

300. 

O.SOGiiOttE 

CO 

-0.827912E 

00 

310. 

0.it21779E 

00 

-0.6U918GF 

00 

320. 

0.352031E 

00 

-0.U700S1E 

00 

330. 

0.296195E 

GO 

-0.2S7025E 

00 

3tt0. 

0  .  2  U  C  9  1 9  E 

00 

-0.839'+';7E- 

-01 

350. 

0. 202 3  3 3E 

00 

0.1372S0E 

00 

125 


HEAVY  WATER-CME    FUEL   ROD 
Z    -    50.    CMS    AXIAL    POSITION 


FREQCCPS) 

AMPLITUDf 

" 

PHASE (RADIANS) 

0. 

0.210  2  9ifE 

'03 

O.OOOOOOE 

00 

10. 

0.20752SE 

03 

0  .  3  2  9  7  8  8  E 

00 

20. 

0.199727E 

03 

0.6531i*3E 

00 

30. 

0.18C175E 

03 

0, 965  0 78 E 

00 

kO. 

0.17ii396E 

03 

0.12631iiE 

01 

50. 

0.159662E 

03 

0.15iiG68E 

01 

60. 

0.1tiH930E 

03 

-0.132527E 

01 

70. 

0.1307li3E 

03 

-0.10C317r 

01 

80. 

0.117'»10E 

03 

-0.822559E 

00 

90. 

0  .  1 0  5  0  3  tt  E 

03 

-0.586925F 

00 

100. 

U.937302E 

02 

-0.360;i61E 

00 

110. 

0.835136E 

02 

-0.1U227f+E 

00 

120. 

0.7it215i*E 

02 

0.GS27StjE- 

-01 

130. 

0.55869iiE 

02 

0  .  2  7 1 6  5  U  E 

09 

I'jO. 

0.533813E 

02 

0.i*63326E 

00 

150. 

0.51692SE 

02 

0.65933t»E 

00 

160. 

0.1*57371E 

02 

0.3'i573iiE 

00 

170. 

0.U0't229E 

02 

0.102693E 

01 

180. 

0.356751E 

02 

0.12039nE 

01 

190. 

0.3U372E 

02 

0.137G22E 

01 

200. 

0.276G33E 

02 

0  . 1  5  'o  S  5  E 

01 

210. 

0  .  2  [;  3  3  2  1 E 

02 

-n.m3ti65E 

01 

220. 

0  .  2  1 3  3  2  9  E 

02 

-0.127461E 

01 

230. 

0.187673E 

02 

-0.111855E 

01 

2'jO. 

0.16i*554E 

02 

-n.954521E 

';o 

2  50. 

0.143779E 

02 

-  0  .  3  1  It  .'t  9  5  E 

00 

250. 

0.125G1(5E 

02 

-0.569060E 

0  0 

270. 

0.10959ftE 

02 

-0.52GQ73E 

CO 

230. 

0.955175E 

01 

-  0  .  5  3  5  1 :.  1  E 

00 

290. 

0  .  8  2  3 1 8  3  E 

01 

-0.2tiG920E 

CO 

300. 

0  .  7 1  3 1 1 G  E 

01 

-0.115i;67E 

CO 

310. 

0.621699E 

01 

0.1G723GE- 

-01 

320. 

0.535539E 

01 

0.142795E 

00 

330. 

0.!;G1.'+32E 

01 

0.265057E 

00 

3U0. 

0.39706I+E 

01 

0.332911E 

0  0 

350. 

G.3tt2235E 

ul 

0.500U65E 

0  0 

125 


FREa(CPS) 

0. 

10. 

20. 

30. 

hO. 

50. 

60. 

70. 

80. 

90. 

100. 

110. 

120. 

130. 

140. 

150. 

160. 

170. 

130. 

190. 

200. 

210. 

220. 

230. 

2U0. 

250. 

260. 

270. 

280. 

290. 

300. 

510. 

320. 

330. 

3ii0. 

550. 

HEAVY  WATER-ONE   FUEL  ROD 
Z    =    55.    CMS    AXIAL    PCS  I Ti ON 

AMPLITUi:)E 


0. 

,17G123E 

03 

0. 

,173G50E 

03 

0. 

,165GU3E 

U3 

0, 

,15627SE 

03 

0. 

,Hi3939E 

03 

0. 

,13093GE 

03 

0. 

,117960E 

03 

0. 

.10553t>E 

03 

0. 

,  9ii0  775E 

02 

0. 

,835550E 

02 

0. 

,  7tj0009E 

02 

0. 

,  5  5  'f  0  3  3  E 

02 

0. 

,  5  7  7  0  6  3  E 

02 

0. 

,503315E 

02 

0. 

,i|ti7122E 

02 

0. 

,392822E 

02 

0. 

,3tiU336E 

02 

0. 

,302275E 

02 

0. 

,26t+695E 

02 

0. 

,23125dE 

02 

0. 

,201823E 

02 

0. 

,176019E 

02 

0. 

,153256E 

02 

0. 

,1335U6E 

02 

0. 

,11G225E 

02 

0. 

,100933E 

02 

0. 

,37537ttE 

01 

0. 

,760165E 

01 

0. 

,659130E 

01 

0. 

,5723^2E 

01 

0. 

,lt969t|0E 

01 

0. 

,l|29173E 

01 

0. 

,367940E 

01 

0. 

.  3 1 3  8 1 5  E 

01 

0. 

,26G15J+E 

01 

0. 

,  2  2  U  2  3  0  E 

01 

PHASE (RADIANS) 

0. 

,OOOOO0E 

00 

0. 

,5i48163E 

00 

0. 

,689599E 

00 

0. 

,101856E 

01 

0. 

,133271E 

01 

-0, 

,1510'i6E 

01 

-0. 

,122721E 

01 

-0. 

,957719E 

00 

-0. 

,700690E 

00 

-0. 

,  t*  5  'f  3  5  9  E 

00 

-0. 

,217o93E 

00 

0. 

,101986E- 

-01 

0. 

,  2  30tf2  5E 

0  0 

0. 

,Uf^3309E 

00 

0. 

,5i»977f|E 

00 

0. 

,8U9810E 

00 

0. 

,10ti471E 

01 

0. 

,  1 2  3  If  6 1|  E 

01 

0. 

,1U1976E 

01 

-0. 

,  1  5 1^  0  9  0  E 

01 

-0. 

,  1 3  G  5  9  0  E 

01 

-0. 

,1191+UOE 

01 

-0. 

,102S1UE 

01 

-0. 

,  3  6  5  6  7 1  E 

00 

-0. 

,705525E 

0  0 

-0, 

,5ti3U01E 

00 

-0. 

,395998E 

00 

-0. 

,2ti729UE 

00 

-0. 

,993772E- 

-91 

0. 

,tiU3227E- 

-01 

0. 

,191S57E 

00 

0. 

,3U1043E 

00 

0. 

,  ti  8  9  3  9  6  E 

00 

0. 

,633517E 

00 

0. 

,775927E 

00 

0. 

,  3 1  3  3  3  0  E 

0  0 
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HEAVY   WATER-ONE   FUEL   ROD 
Z    =    60.    CMS    AXIAL    POSITION 


FREQCCPS) 

AMPLITUDE 

PliASE(RADIAMS) 

0. 

0.14637[iE 

03 

O.OOOOOOE  00 

10. 

O.lUiilGSE 

03 

0.367907E  00 

20. 

0.1379S1E 

03 

0.72S379E  00 

30. 

0.123878E 

05 

0.10759t;E  01 

hO. 

0.1180G6E 

03 

O.lhOlkZE    01 

50. 

0.106693E 

03 

-O.U1952E  01 

50. 

0  .  9  5  U  U  3  7  F 

02 

-0.112117E  01 

70. 

0.Sli7998E 

02 

-0.837SStjE  CO 

80. 

0.7l»985SF 

02 

-0.568002E  00 

90. 

0.560365E 

02 

-0. 31007 '^E  00 

100. 

0  .  5  8 1 1 2  1 E 

02 

-U.625ti68E-01 

110. 

0.50992I1E 

02 

0.175375F  00 

120. 

0.iiUG719E 

02 

0.ii05372E  00 

130. 

0.390710E 

02 

0.623  It  SOE  00 

lUO. 

0.3ttl210E 

02 

0.8ii3692E  00 

150. 

0.297752E 

02 

0.105253E  01 

150. 

0.259U72F 

02 

0. 12 5 5 75 F  01 

170. 

0.225777E 

0  2 

0.1U5330r  01 

180. 

0.196221E 

02 

-0.1U365I+E  01 

190. 

0.170U31E 

02 

-0.13 096 3E  01 

200. 

0.1I+7779E 

02 

-0.112759E  01 

210. 

0.12S138E 

02 

-0.95037 IE  00 

220. 

0.111013E 

02 

-0.777295E  CO 

230. 

0.961035E 

01 

-0.60356nE  00 

2I4O. 

0.831SGCE 

01 

-0.41*168  IE  0  0 

250. 

0.713039E 

01 

-0.2795U6E  00 

260. 

0. 6 2 1 5 0 2  F 

01 

-0.120307E  00 

2  70. 

0.536257F 

01 

0.399535E-01 

230. 

0.fiG17UlE 

01 

U.195G10F  0  0 

290. 

0.397693E 

01 

0.350521E  00 

300. 

0  .  3  J*  2  '^  U  9  E 

01 

0.50G518E  00 

510. 

0.293166E 

01 

0.669156E  00 

520. 

0.2i47711F 

01 

0.32S590F  00 

330. 

0.20S31GE 

01 

0.933026E  00 

3I1O. 

0.173539F 

01 

0.11 3315 F  01 

3  50. 

0.1t+U232E 

01 

0.127171E  01 
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HEAVY  WATER-ONE    FUEL   ROD 
Z    =    55.    CI'iS    AXIAL    POSITIOI! 


FREQ(CPS)  AMPLITUDE 

0.  0.117520E  03 

10.  0.1156't8E  03 

20.  0.110U27E  03 

30.  0.102757E  03 

hO.  0.9372  9  7E02 

50.  0.8U2U5ltE  02 

60.  0.749378n  02 

70.  0.66180UE  02 

80.  0.581679r:  02 

00.  0.509578E  02 

100.  0.l|ii5375E  02 

110.       •  0.388603E  02 

120.  0.333650E  02 

130.  0.29ii773E02 

lliO.  0.256277E  02 

150.  0.222G16E  02 

160.  0.193173E  02 

170.  0.1573U7E  02 

180.  0.1'|i4S36E  02 

190.  G.12520itE  02 

20  0.  0.10C297E02 

210.  0.937138E  01 

220.  0.810917E01 

230.  0.700!+15E  01 

2k0.  0.503ii52E  01 

250.  0.516756E  01 

260.  0.t»'t2165E  01 

270.  0.378505E  01 

230.  0.321975E  01 

290.  0.273707E  01 

300.  0.230132E  01 

310.  0.19U3nSE  01 

320.  0.163  9  9  5E01 

330.  0.1395 8 BE  01 

3U0.  0.121313E  01 

35  0.  0.105133E  01 
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HEAVY   WATER-ONE   FUEL   ROD 
Z    =    70.    CMS    AXIAL    POSITION 


FRF.Q(CPS) 

AMPLITUDE 

PHASE (RADIANS) 

0. 

0.989173E 

02 

O.OOOOOOE 

00 

10. 

0.972562E 

02 

0.1+0721+2E 

00 

20. 

0.926077E 

02 

0.806277E 

00 

30. 

0.35799SE 

02 

0.119086E 

01 

40. 

0.773036E 

02 

0.155736E 

01 

50. 

0.69l»506E 

02 

..     -0.123671E 

01 

60. 

0.613053E 

0  2 

-0.907707E 

00 

70. 

0.537035E 

02 

-0.595735E 

00 

80. 

0.1+0  8  050E 

02 

-0.29G978E 

00 

90. 

0.t+063S6E 

02 

-0.15S363E- 

■01 

100. 

0.551921+E 

02 

0  .  2  5  5  2  2  1 E 

00 

110. 

0.301+16  7E 

02 

0.515U86E 

00 

120. 

0.262508E 

02 

0.765955E 

00 

130. 

0.225307E 

02 

0.100797E 

01 

ItiO. 

0.19it77GE 

02 

0  . 1 2  1+  2  0  5  E 

01 

150. 

0.167I+71E 

02 

0  .  1 1*  5  8 1+ 1+  E 

01 

160. 

0.11+3735E 

02 

-  0  . 1 1+  5  3  9  0  E 

01 

170. 

0.123'+7GE 

02 

-0.121+21+7E 

01 

ISO. 

O.IOGIOI+E 

02 

-0.103523E 

01 

100. 

0.9 0992 3E 

01 

-0.832S50E 

00 

200. 

0. 779551E 

01 

-0.635099E 

00 

210. 

0.66G195E 

01 

-  0  . 1+ 1+  2  5 1+  3  £ 

00 

220. 

0.569325E 

01 

-0.25.'+790E 

oc 

230. 

0.1+SU521E 

01 

-0.7001+90E- 

-01 

21+0. 

0.I+11175E 

01 

0.10G9I+1E 

00 

250. 

0.3it8390E 

01 

U.277253E 

00 

260. 

0.295708E 

01 

O.U39ig8E 

00 

270. 

0.251I+51E 

01 

0  .  5  9  9  0 1+  9  E 

00 

2  30. 

0  .  2  1 3  7  0  6  E 

01 

0.752  515E 

00 

290. 

0.1821+91+E 

01 

0.9051+18E 

00 

300. 

0.15i+7ijlE 

01 

0.105S18E 

01 

310. 

0.131233E 

01 

0.119631+E 

01 

320. 

0  .  11 2  9  3  0  E 

01 

0.1333UE 

01 

330. 

0.97U937E 

00 

0.i:+7i+99E 

01 

31+0. 

0.853550E 

00 

-0.151556E 

01 

350. 

0. 737710E 

00 

-9.133359E 

ni 
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HEAVY  WATER-ONE   FUEL   ROD 
I    =    75.     CMS    AXIAL    POSITION 


FREQ(CPS) 

AMPLITUDE 

PIIASE(RADIAr;S) 

0. 

0.820f*';lE    02 

O.OOOOOOE    00 

10. 

0.806118E    02 

0.U25185E    00 

20. 

0.7659U3E    02 

0.8ii2028E    00 

30. 

0.707103E    02 

0.12i*378E    01 

kO. 

0.633256E    02 

-0.151'»7itE    01 

50. 

0.555if63E    02 

-0.115151E    01 

60. 

0.U90807E    02 

-0.807876E    00 

70. 

0.t»32239E    02 

-0.tt82102E    00 

80. 

0.373900E    02 

-0.172515E    00 

90. 

0.322265E    02 

0.122'+92E    00 

100. 

0.277033E    02 

0.lt05058E    00 

110. 

0.237605E    02 

0.676320E    00 

120. 

0.203'|it3E    02 

0.937229E    00 

130. 

0.173950E    02 

0.118382E    01 

lUO. 

0.U3532E    02 

0.1i»3176E    01 

150. 

0.120639E    02 

-0.1£t7520E    01 

luO. 

0.107937E    02 

-0.12'48[t5E    01 

170. 

0.920699E    01 

-0.10239t4E    01 

130. 

0. 735106E    01 

-0.81U653E    00 

19  0. 

0.663751E    01 

-0.e053i4^E    00 

200. 

0.563723E    01 

-0.I40072GE    00 

210. 

0.tiS25  0  2E    01 

-0.201390E    00 

220. 

0.I40C763E    01 

-0.923031E-02 

230. 

0.3'4G150E    01 

0.176200E    0  0 

2U0. 

0.293757E    01 

0.557092E    00 

250. 

0.2ti93U7E    01 

0.53U927E    00 

260. 

0.211685E    01 

0.712595E    00 

270. 

0.173i;98E    01 

0.886961E    00 

280. 

0.15077I4E    01 

0.105721E    01 

290. 

0.125739E    01 

0. 122659 E    01 

300. 

0.10506GE    01 

0.137';73E    01 

310. 

0.883S12E    00 

0.1525filE    01 

320. 

0.7't572i+E    00 

-O.Hi7559E    01 

330. 

0.633579E    00 

-0.132U15E    01 

jliO. 

0.53't991E    00 

-O.llSOO'+E    01 

3  50. 

0,!;55317E    00 

-0. ] niS3nE    CI 
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HEAVY   WATER-ONE   FUEL   ROD 
80.    ens    AXIAL    POSITION 


FREQ(CPS) 

AMPLITUDE 

PHASE(RADIANS 

0. 

0.677836E 

02 

O.OOOOOOE 

00 

10. 

0.665593E 

02 

O.U[tii392E 

00 

20. 

0.531157E 

02 

0.8802(j3E 

00 

30. 

0.580734E 

02 

0.13095GE 

01 

40. 

0  .  5  2  1 8  8  (4  E 

02 

-0.1tt(»017E 

01 

50. 

0.iiG07';5E 

02 

.  -0.106015E 

01 

60. 

0.1*01696E 

02 

-0.700599E 

CO 

70. 

0.3i>7170E 

02 

-0.359893E 

00 

8U. 

0.293315E 

02 

-0.3G5220E- 

-01 

00. 

0. 255322E 

02 

0.271527E 

0  0 

100. 

0.217905E 

0.565950E 

0  0 

lie. 

0.18553GE 

02 

0.8tt73.'*3E 

GO 

12  0. 

0.1578S7E 

02 

0.111850E 

01 

130. 

0.13'a79E 

02 

0.137399E 

01 

liiO. 

0.111>013E 

02 

-0.151170E 

01 

150. 

0.963S70E 

01 

-0.12G9Ur 

01 

IGO. 

0.32Ult|6E 

01 

-0.103367E 

01 

170. 

C.700139E 

01 

-0.803113E 

00 

180. 

0.59iill6E 

01 

-0.5S0755E 

00 

190. 

0.50U773E 

01 

-0.3G2635E 

00 

200. 

0.it275GGE 

01 

-0.150117E 

00 

210. 

0.362053E 

01 

0.557071E- 

■01 

220. 

0.30G131E 

01 

0.2550S1E 

00 

230. 

0.259400E 

01 

0.lt50272E 

00 

2U0. 

0.2192I42E 

01 

O.ettlUOOE 

oc 

250. 

0.185763E 

01 

0.82!i'i23E 

0  0 

2G0. 

0.157258E 

01 

0.100817F 

01 

270. 

0.133052E 

01 

0.11S205E 

01 

280. 

0.113077E 

01 

0.135897E 

01 

290. 

0.952901E 

00 

0.15317SE 

01 

300. 

G.801911E 

00 

-0,lkhhhh^. 

01 

310. 

0.67G259E 

00 

-0.129336E 

01 

320. 

0.577509E 

00 

-0.i:3917E 

01 

330. 

O.US823iiE 

00 

-0.993791E 

00 

3'*0. 

0.[t2'tG73E 

00 

-0.86(i031E 

00 

350. 

0.5  77957E 

00 

-0.799!4l2E 

0  0 
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HEAVY  WATER-ONE   FUEL   ROD 
Z    =    85.    CfIS    AXIAL    POSITION 


FRCQCCPS) 

AMPLITUDE 

PHASE  C^.ADIAriS) 

0. 

0.5G059CE 

02 

O.OOOOOOE 

00 

10. 

0.55002  i»E 

02 

0.l;bl676E 

09 

20. 

0.520537E 

02 

0.914350E 

00 

30. 

0.I»77564E 

02 

0.135110E 

01 

'jO. 

0.it27fi'tttE 

02 

-0.157557E 

01 

50. 

0.3755'*3E 

02 

-0.973156E 

00 

GO. 

0.32552fiE 

02 

-0.603752E 

00 

70. 

0.27956SE 

02 

-0.2ti9117E 

00 

80. 

0.233516E 

02 

0.875't07E- 

-01 

90. 

0.20257i»E 

02 

0.lt07696E 

00 

100. 

0.1716lt3E 

02 

0.71320tiE 

00 

110. 

0.1f*5120E 

02 

0.100500E 

01 

120. 

0.122575E 

02 

0.12 355 7E 

01 

130. 

O.IO3I471+E 

02 

0.15573UE 

01 

11*0. 

0.8725147E 

01 

-0.13227UE 

01 

150. 

0.7350(j5E 

01 

-0.1069G3E 

01 

160. 

0.613387E 

01 

-0.825Sti0E 

00 

170. 

0.5203ii8E 

01 

-0.531030E 

00 

180. 

0.U38G51E 

01 

-0.363976E 

00 

190. 

0.370'f26E 

01 

-0.1141777E 

00 

200. 

0.312598E 

01 

0.757S33E- 

-01 

210. 

0  .  2  6  3  8  7  G  E 

01 

0.283070E 

00 

220. 

0.2226g6E 

01 

0. 1*96999  E 

00 

230. 

0.1875'+3E 

01 

0.701503E 

0  0 

2U0. 

0  .  1 5  3 1 '»  7  E 

01 

0.899952E 

00 

250. 

0.133523E 

01 

0  . 1 1 9 1  ?  2  E 

01 

260. 

0.111615E 

01 

0.130267E 

01 

270. 

0.928116E 

00 

0.1U9395E 

01 

280. 

0.767007E 

00 

-0.1t;6015E 

01 

290. 

0.65161UE 

ou 

-0.123900E 

01 

300. 

0.525G77E 

00 

-0.112299E 

01 

310. 

0.!i3520SE 

00 

-0.955it95E 

CO 

320. 

0.361311E 

00 

-0.793715E 

00 

330. 

0.299702E 

00 

-0.51S755E 

00 

3ft0. 

0.2ft0t{95E 

00 

-0.ft33855E 

90 

350. 

0.1S5072E 

00 

-0.279697E 

CO 

APPEiroiX   F 
NU'^RICAL  CONSTAIiTS 
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The  natural  uranium  fuel  element  used  in  this  experiment  has  been 
described  in  detail  by  Dunlap  (6).  The  fuel  rod  consisted  of  six  cylin- 
drical slugn  contained  in  an  aluminum  tube.   Numerical  values  used  for 
other  constants  were  as  follows . 

Thermal  Diffusion  Constant  (D^,  ) O.815I4  cm 

th 

Diffusion  Length  Squared  (^  +v^) 

2 
99'0'A   Heavy  Water 3738  era 

99.5^  Heavy  Water 6751  cm^ 

Thermal  Velocity  ( v  ^  ) 2.k^   x  10''  cm/sec 

th 

Age  to  Thermal  (f.,  ) 

99.0^  Heavy  Water II8.8  cm 

2 
99 . 5'/i  Heavy  Water 119 .  ^   cm 

Slowing  Down  Time  (L  ) 5.08  x  10~  sec 

Resonance  Escape  Probability  (p) .999 

Fast  Fission  Factor  (t^  ) 2.07 

2  _3   2 

Fundaiiiental  Buckling  (3   ) I.263  x  10   cm 

First  Order  Buckling  (B^) 6. 6^2  x  lO"^  cm 

Extrapolated  Radius  of  System  (?•„„„) 67  .67  en 

Radius  of  Fuel  Rod  (r  ) 1.27  cm 

o 


1.   Heavy  water  parameters  were  obtained  from  II.    J.  Diaz  by  private 
communication. 
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